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ABSTRACT 


HYDROGRAPH  RISE  TIMES 

The  runoff  distribution  with  respect  to  time  is  needed  for  the 
design  of  many  hydraulic  structures.  In  most  of  the  presently  avail¬ 
able  techniques  of  flood  peak  estimation  and  design  hydrograph  pre¬ 
diction,  it  is  necessary  to  determine  characteristic  or  critical 
hydrograph  times.  This  study  suggests  the  use  of  the  hydrograph 
rise  time,  which  is  the  time  between  the  commencement  of  runoff 
and  the  maximum  discharge  measured  at  a  stream  gaging  station, 
as  the  time  characteristic  of  a  particular  watershed.  It  is  shown, 
however,  that  rise  times  within  a  watershed  have  considerable 
variability  but  a  general  distribution  of  values  may  be  applied  to  all 
the  watersheds  used  in  this  study.  Because  of  the  inadequacy  of  the 
existing  rise  time  prediction  methods,  several  new  equations  are 
derived  in  a  stepwise  regression  analysis  by  computer.  The  inves¬ 
tigation  uses  407  flood  events  of  47  watersheds  in  1  3  states,  which 
is  a  larger  sample  than  the  data  in  any  previous  study  of  the  same 
type.  The  conclusions  are  limited  to  floods  caused  by  thunderstorms 
or  storms  of  short  duration  (two  or  three  hours)  which  are  generally 
the  most  critical  for  small  watersheds.  Different  equations  are 


in 
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obtained  for  estimating  median  rise  times  in  humid  regions  and  arid 
regions. 
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RESUME 


Les  Temps  de  Montee 

Dans  la  conception  de  nombreuses  structures  hydrauliques 

on  a  besoin  de  connaitre  la  repartition  de  l'ecoulement  dans  le  temps. 

* 

Les  techniques  existantes  pour  la  determination  de  la  pointe  de  crue 
et  pour  la  prediction  de  l'hydrogramme  requierent  le  calcul  des 
temps  caracteristiques  ou  critiques  de  l'hydrogramme.  On  a  suggere 
l'usage  du  temps  de  montee  comme  le  temps  caracteristique  d'un 
bassin  versant.  Ce  temps  represente  la  periode  allant  du  commen- 
cement  de  l'ecoulement  jusqu'a  la  pointe  de  crue.  Cependant,  on  a 
montre  que  les  temps  de  montee  dans  un  bassin  varient  considerable - 

ment,  mais,  qu'une  distribution  generate  de  valeurs  peut  etre 

-  *•  ^  ^  -• 
appliquee  a  tous  les  bassins  etudies.  Comme  les  methodes  existantes 

■-  >. 

pour  determiner  les  temps  de  montee  ne  sont  pas  adequates  on  a 

>*  ^  >* 
derive  plusieurs  formules  nouvelles  par  analyse  regressive  en 

> 

employant  un  calculateur  electronique.  Dans  cet  etude  on  a  utilise 

407  crues  dans  47  bassins  et  13  etats  -  c'est  le  plus  grand  echantillon 
-*  *  * 

de  donnees  utilise  dans  des  etudes  de  ce  genre.  Les  conclusions  sont 
limitees  aux  crues  causees  par  les  averses  ou  les  pluies  de  courte 
duree  (deux  ou  trois  heures  )  qui  sont  generalement  les  plus  critiques 
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dans  les  petits  bassins.  On  a  obtenu  des  formules  differentes  pour 

>  * 

calculer  la  mediane  des  temps  de  montee  dans  les  regions  humides 
et  dans  les  regions  arides. 
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Chapter  I 


INTRODUCTION 

Numerous  methods  are  presently  available  to  estimate  the 
total  volume  of  runoff  associated  with  a  particular  combination  of 
storm  rainfall  and  antecedent  conditions  (1,  2,  3)  *  .  In  many  in¬ 
stances,  however,  the  computation  of  runoff  volume  is  only  a 
preliminary  step  in  the  determination  of  the  outflow  hydrograph 
from  a  basin.  The  design  of  hydraulic  structures  may  be  cited  as 
a  case  where  the  engineer  is  concerned  not  only  with  the  total  volume 
but  also  with  the  maximum  discharge  and  its  time  of  occurrence. 

Various  studies  have  provided  methods  of  estimating  flood 
peaks  directly  from  watershed  characteristics  such  as  area,  slope, 
location  etc.  .  .  (4,  5,  6).  It  appears  possible  to  derive  similar  rela¬ 
tionships  for  estimating  rise  time  and  recession  characteristics 
so  that  the  complete  design  hydrograph  may  be  synthesized. 

For  small  watersheds,  flood  frequencies  are  often  derived 
from  rainfall  frequencies  by  techniques  such  as  Unit-graph  method 
(7),  "Rational"  Formula  (8),  Soil  Conservation  Service  Hydrograph 
Families  (9),  Bureau  of  Public  Roads  method  (10),  Chow's  method 

*  Numbers  in  the  parentheses  refer  to  the  references  beginning  on 
page  80. 
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(11),  and  Tacitly  Maximized  Peaks  (12).  In  most  of  these  it  is  ne¬ 
cessary  to  estimate  characteristic  or  critical  hydrograph  times  that 
are  related  to  the  rise  time  and  effective  rainfall  durations.  These 
are  mainly  watershed  characteristics,  long  characteristic  times 
being  associated  with  large  or  slow-responding  watersheds  and  vice 
versa. 


Possibly  no  part  of  the  hydrologic  literature  is  more  con¬ 
fusing  than  that  devoted  to  the  determination  of  these  times.  There 
are  many  slightly  different  concepts,  some  of  which  have  conflicting 


interpretations.  They  include  "rise  time",  "storage  delay",  "lag". 


M  1 1 


II  lU _ II 


"critical  duration",  "effective  duration",  "optimum  duration"  and 


it. 


time  to  peak".  Further,  there  are  many  suggested  methods  of 


estimating  these  values  but  most,  if  not  all,  are  of  uncertain 
reliability  (13).  This  study  should  clarify  the  significance  of  the 


concepts  and  provide  some  recommendations  for  improved  estimates 


of  the  required  times. 


. 

» 
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Chapter  II 

REVIEW  OF  LITERATURE 

As  mentioned  in  the  introduction,  no  part  of  the  hydrologic 
literature  is  more  confusing  than  that  devoted  to  the  determination 
of  time  parameters.  This  chapter  is  limited  only  to  the  literature 
survey  on  rise  time  or  its  equivalent. 

1 .  Ramser,  1927  (14) 

Ramser  conducted  a  series  of  rainfall  and  runoff  measure¬ 
ments  on  six  agricultural  watersheds  ranging  in  area  from  1.  25 
acres  to  112.0  acres,  near  Jackson,  Tennessee.  The  time  of  con¬ 
centration  was  defined  as  the  time  required  for  the  water  to  flow 
from  the  farthest  point  on  the  watershed  to  the  gaging  station.  This 
was  determined  by  noting  the  rise  time  which  is  the  time  required 
for  the  water  in  the  channel  at  the  gaging  station  to  rise  from  the 
low  to  the  maximum  stage  as  recorded  by  the  water-stage  recorders. 
Ramser  found  that  this  period  varied  to  some  extent  for  the  different 
rains,  depending  upon  the  degree  of  saturation  of  the  watersheds  at 
the  occurrence  of  the  rain  that  produced  the  maximum  rate  of  runoff. 
He  also  found  that  the  time  of  concentration  would  be  less  if  the 


channel  were  partially  filled  with  water  when  the  rains  of  greatest 


. 
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intensity  occurred,  and  would  be  less  for  rains  of  high  intensity  than 
for  rains  of  less  intensity.  However,  owing  to  the  inability  of  the 
recording  instruments  to  register  the  time  with  high  accuracy,  the 
time  of  concentration  or  rise  time  was  taken  to  be  the  same  for  the 
different  rains  on  the  same  watershed. 


2.  Kirpich,  1940  (20) 

Kirpich  presented  two  curves  based  on  data  obtained  by 
Ramser  (14)  for  evaluation  of  time  of  concentration.  In  the  two 
curves  shown  (figures  1  and  2)  it  is  assumed  that  the  time  of  con¬ 
centration  for  any  watershed  depends  on  a  factor,  K,  which  varies 
directly  with  the  length  of  travel  and  inversely,  with  the  square 
root  of  the  slope.  It  could  be  put  into  a  form  of  empirical  formula 
(21  )  as  follows: 


T 

C 


0. 00013 


L 

s 


0.  77 
0.  385 


where  L  is  the  length  of  the  basin  areas  in  miles,  measured  along 
the  water  course  from  the  gaging  station  and  in  a  direct  line  from 
the  upper  end  of  the  watercourse  to  the  farthest  point  on  the  drainage 
basin;  and  S  is  the  ratio  in  feet  to  L  of  the  fall  of  the  basin  from  the 
farthest  point  on  the  basin  to  the  outlet  of  runoff,  or  approximately 
the  average  slope  of  the  basin  in  dimensionless  ratio. 

Considering  the  original  data  obtained  by  Ramser  (14),  this 
T  is  actually  the  rise  time. 


. 

> 

* 
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Values  of  K 


D 


Figure  1.  Concentration  time 
of  small  agricultural  water¬ 
sheds,  based  on  K 
(  after  Z.  P.  Kirpich,  1940) 
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Figure  2.  Concentration  time 
of  small  agricultural  water¬ 
sheds,  based  on  K 
(  after  Z.P .  Kirpich  ,  1940) 


is  the  difference  in  elevation  in  feet  between  the  most  remote 
point  and  the  outlet. 

H  is  the  average  height  of  the  watershed  above  the  outlet  in  feet. 
L  is  the  length  of  travel  in  feet. 

A  is  the  area  of  the  watershed  in  square  miles. 


These  definitions  are  referred  to  the  work  of  Kirpich  only. 


6 


3.  Soil  Conservation  Service,  1957  (9) 

The  Soil  Conservation  Service  developed  a  nomograph,  figure 
3,  from  Kirpich's  paper  (20),  which  is  based  on  Ramser's  experi¬ 
mental  data  (14).  In  the  original  study  of  Ramser's,  the  values  pre¬ 
sented  were  actually  the  hydrograph  rise  times  although  they  were 
called  times  of  concentration.  The  SCS  has  apparently  overlooked 
this  because  in  their  report  they  regard  the  times  to  peak  as  different 
to  Ramser's  original  times  of  concentration.  The  SCS  nomograph  is 
suggested  for  estimating  the  time  of  concentration  T^,  in  hours  knowing 
the  length  in  feet,  L,  of  the  longest  waterway  from  the  watershed  out¬ 
let  to  the  ridge  and  the  difference  in  elevation,  H,  between  the 
watershed  outlet  and  the  farthermost  point,  in  feet. 

The  SCS  suggested  that  an  average  time  of  concentration  can 

be  obtained  using  the  time  from  beginning  of  rise  to  peak,  T  ,  in  the 

R 

equation 

tr=  tc°'5+0'6Tc 

Figure  4  permits  an  easy  solution  for  T^  once  the  T^  is  known.  T^ 
is  averaged  value  for  simple  hydrographs  or  from  field  observations 
of  simple  hydrographs. 

For  the  design  hydrograph,  the  SCS  suggested  to  estimate  the 

time  to  peak,  T  ,  by  relating  it  to  the  duration  of  rainfall  excess,  D  , 

and  the  time  of  concentration,  T^,  of  the  watershed. 

T  =  0.  5  D  +0.6  T 
R  e  C 


; 

k'  >  ri  A  i »f  r*l  , 
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L  =  7,250  feet 
H  =  130  feet 
then  Tc  =0.57  hours 
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Figure  3.  Time  of  concentration  (after  SCS-USDA,  1957) 
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Figure  4.  Time  to  peak  (after  SCS-USDA,  1957) 
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4.  Minshall,  1960  (1  ) 

Using  three  Edwardsville  watersheds  in  Illinois  with  a  large 
number  of  runoff  periods  each  year,  Minshall  presented  a  method 
for  constructing  synthetic  unit  hydrographs  for  drainage  areas  from 
20  acres  up  to  about  500  acres  and  for  different  rainfall  intensities. 

Time  from  beginning  of  rainfall  excess  to  peak  rate  of  runoff, 
Tp,,  can  be  obtained  from  the  nomograph  in  figure  5  for  a  particular 
size  area  and  different  rainfall  intensities. 


5.  Gray,  1961  (15) 

In  1961,  Gray  found  in  his  study  of  lag  time  for  drainage  areas, 
ranging  in  size  from  25  to  50  square  miles  in  the  central  United  States, 
that  time  lag  could  be  correlated  to  the  time  period  of  rise  of  the  hy¬ 
drograph.  From  94  selected  storms  he  concluded  that  : 


T  =  0.  996  T 
-L  R 


1. 005 


where  was  defined  as  lag  time  from  center  of  mass  of  rainfall  to 
peak  of  runoff,  in  hours,  and  T  was  the  time  period  of  rise  of  the 
hydrograph.  He  stated  that  since  the  coefficient  and  exponent  were 
sufficiently  close  to  unity,  they  could  be  considered  numerically 


equal  to  one  and  thus  "a  given  change  in  T  produces  an  equal  change 
in  T  ".  Since  the  time  of  rise  was  found  to  be  equal  to  the  lag  time, 

1j 

he  stated  that  T  could  be  used  as  an  important  variable  in  relating 

R 


the  hyetograph  to  the  hydrograph. 


- 
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Figure  5.  Effect  of  rainfall  intensity  and  size  of  watershed  on  the  unit  hydrograph 

(after  N.  E.  Minshall,  1960) 
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6.  Keppel  and  Renard,  1962  (16  ) 

In  a  paper  by  Keppel  and  Renard  an  equation  for  average  rise 
time  for  the  Walnut  Gulch  watersheds  was  given  as 


25-3  AW 


-0.  14 


where: 


T  =  rise  time  in  minutes 
K 

=  watershed  area  in  square  miles. 

As  explained  in  that  paper,  the  decrease  in  rise  time  with  increase 
of  drainage  area  is  probably  due  to  two  interrelated  factors:  the 
occurrence  of  most  of  the  transmission  losses  during  rising  stages 
of  the  flow;  and  the  presence  of  overriding  translatory  waves  as  the 
flow  moves  through  the  channel. 


7.  Chow,  1962  (1  1  ) 

Chow,  after  analyzing  the  data  from  20  small  midwestern 
watersheds,  recommended  time  to  peak  be  determined  by  the  equa¬ 
tion: 

T  =  0.  00236  (L/S  °*  5)  °'  64 
R 

In  the  foregoing  L  is  the  channel  length  in  feet  as  measured  up  the 
main  water  course  to  the  divide  and  S  is  the  slope  in  per  cent.  The 
latter  is  determined  by  plotting  the  stream  profile  and  fitting  a 
straight  line  through  the  gaging  station  so  that  the  area  between  the 


' 

I 

* 

►  c 


line  and  the  profile  lying  below  the  line  is  equal  to  that  lying  above  it. 

T  is  measured  in  hours,  and  assumed  to  be  constant  for  each  water¬ 
shed. 

8.  Wu,  1963  (17  ) 

I-Pai  Wu  made  a  study  to  analyze  the  existing  data  on  21  small 

watersheds  (less  than  100  square  miles)  distributed  throughout  the 

state  of  Indiana  and  to  determine  the  relations  between  the  slope  of 

the  hydrograph  and  some  identifiable  and  readily  obtainable  watershed 

characteristics.  It  was  noted  that  the  time  to  peak,  T  ,  did  not  vary 

R 

radically  for  the  same  watershed  and  hence  an  average  time  to  peak 
was  used. 

He  found  that  average  time  to  peak,  T  ,  is  correlated  with 
three  measurable  watershed  characteristics,  A  (watershed  area), 

L  (length  of  main  stream),  and  S  (mean  slope  of  main  stream). 

T  =  31.42  A1'085  L/1'233  S'0'668 
R 

Wu  suggested  that  the  design  storm  duration  be  made  equal  to 

the  value  of  T  .  The  selection  of  the  duration  of  storm  rainfall  is 

K 

based  on  the  time  of  concentration  of  the  watershed  considering  that  a 
storm  rainfall  duration  equal  to  the  time  of  concentration  will  result  in 
the  maximum  rate  of  discharge.  Since  the  time  of  concentration  of  the 
watershed  is  difficult  to  determine,  it  may  be  assumed  that  the  design 
rainfall  is  uniformly  distributed  over  the  entire  watershed,  in  which 


r  •  l  ■■  9  1  , 
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case  the  time  of  concentration  is  equal  to  the  time  to  peak. 


9.  Kibler,  1965  (18  ) 

Using  a  total  of  140  selected  storms  pertaining  to  runoff  from 
the  Shaver  Creek  basin,  a  3.7  5  square  mile  watershed  in  central 
Pennsylvania,  Kibler  attempted  an  emperical  approach.  By  the  me¬ 
thod  of  multiple  regression,  he  developed  the  following  equation: 


T 


R 


5. 538  D 


0.  333 


in  which  T  is  defined  as  total  time  elapsed  between  initiation  of  sur- 
R 

face  runoff  and  time  of  peak  discharge,  in  hours,  and  D  is  the  storm 
duration  in  hours.  This  equation  is  considered  valid  for  prediction 
of  time  to  peak  of  the  March- May  storms  with  average  errors  of  17.9 
per  cent.  It  is,  however,  limited  to  the  Shaver  Creek  basin. 


1 0.  Machmeier,  1966  (19  ) 

Machmeier  developed  a  mathematical  model  of  a  .21.  35  square 
mile  watershed  with  the  land  and  channel  characteristics  representa¬ 
tive  of  small  watersheds  in  southeastern  Minnesota.  In  one  part  of 
his  study  he  attempted  to  show  how  the  various  time  parameters  are 
affected  by  supply  rate  and  duration. 

For  a  duration  of  one  hour,  figure  6  shows  the  effect  of  the 
supply  rate  on  the  time  to  peak.  For  each  of  the  durations  used,  the 
time  to  peak  for  the  model  decreases  as  the  supply  rate  increases. 


f  .!  '  1  V.  " 
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Figure  6.  Model  watershed  hydrographs  for  various  supply  rates 

and  a  duration  of  1 .  0  hour  (after  R.  E.  Machmeier,  1966) 
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Figure  7.  Model  watershed  hydrographs  for  various  durations  and 

a  supply  rate  of  1 .  0  in. /hr  (after  R.  E.  Machmeier,  1966) 
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Machmeier  explained  that  this  nonlinear  response  is  evidently  due  to 
the  fact  that  the  average  velocity  in  the  channel  system  increases  with 
discharge.  According  to  unit  hydrograph  theory,  the  time  to  peak  is 
not  affected  by  supply  rate. 

The  effect  of  duration  on  time  to  peak  is  shown  in  figure  7  for 
a  supply  rate  of  1.  00  inch  per  hour.  Generally,  one  would  expect 
time  to  peak  to  increase  with  duration.  However,  the  results  of  this 
study  indicate  that  time  to  peak  has  a  minimum  value  and  tends  to 
increase  both  for  longer  and  shorter  durations.  This  same  trend  was 
exhibited  for  all  of  the  supply  rates  tested. 

All  of  the  hydrologic  type  equations  must  remain  approxima¬ 
tions  and  individual  watersheds  will  depart  appreciably  from  mean 
relationships  presented  previously.  Among  the  factors  not  easily 
evaluated  in  analyzing  existing  data  are:  the  relative  role  of  overland 
and  channel  flow,  the  effect  of  intensity  and  duration  of  rainfall,  the 
effect  of  roughness  on  both  overland  and  channel  flow,  and  the  var¬ 
iance  in  basin  morphology  as  represented  by  such  items  as  drainage 
density,  channel  and  land  slopes,  and  channel  cross-section  through¬ 


out  the  watershed. 
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CHAPTER  III 

ANALYSIS  AND  INTERPRETATION  OF 
HYDROGRAPH  RISE  TIME 

Inability  to  measure  the  time  most  desired  by  theoretical 
definition  is  one  of  many  reasons  for  confusion  in  the  determination 
of  characteristic  times.  Rise  time,  sometimes  referred  to  as  time 
to  peak,  marks  the  occurrence  of  the  runoff  peak  rate;  hence,  it  has 
a  definite  physical  significance  and  can  be  measured  easily.  There¬ 
fore,  the  rise  time  will  be  used  in  this  study  as  a  main-time  para¬ 
meter  (or  characteristic).  When  necessary,  it  can  be  related  to 
other  hydrograph  characteristics  for  predicting  the  time  distribution 
of  runoff. 

1 .  Hydrograph  Characteristics 

The  hydrograph  can  be  regarded  as  an  integral  expression  of 
the  physiographic  and  climatic  characteristics  that  govern  the  rela¬ 
tions  between  rainfall  and  runoff  of  a  particular  drainage  basin  (21). 

It  shows  the  time  distribution  of  runoff  at  the  point  of  measurement, 
defining  the  combined  effects  of  the  rainfall  excess  and  the  complexities 
of  the  basin  characteristics  by  a  simple  empirical  curve. 

A  typical  single-peaked  hydrograph,  figure  8,  consists  of  three 
parts:  the  approach  segment  AB,  the  rising  segment  BC  and  the  falling 
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Figure  8.  A  Typical  hydrograph  with  corresponding  hyetograph 
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Figure  9.  Definition  sketch  of  selected  measures  of  time  scale 
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segment  CD.  The  approach  segment  represents  antecedent  flow  before 
significant  runoff-producing  rainfall  reaches  the  outlet.  The  rising 
segment  usually  commences  immediately  after  the  runoff-producing 
rainfall  commences.  It  represents  increasing  rates  of  runoff  caused 
by  the  inflow  to  channel  storages  from  the  rainfall.  The  peak  of  a 
hydrograph  is  the  highest  concentration  of  the  runoff  from  a  drainage 
basin.  The  falling  segment  represents  withdrawal  of  water  from 
storage  after  the  main  rainfall  excess  has  finished.  It  is  usually 
considered  to  be  largely  independent  of  the  time  pattern  of  rainfall 
excess. 

2.  Relationship  Between  Rise  Time  and  Other  Hydrograph 

Characteristics 

The  relationships  between  rise  time  and  other  hydrograph 

characteristics  are  illustrated  in  figure  9  where  T^  is  the  time 

difference  between  the  lowest  rate  of  runoff  q.  and  the  peak  rate  q 

l  o 

of  the  hydrograph.  It  is  assumed  that  the  rainfall  excess  is  of 
uniform  intensity  and  its  time  of  commencement  corresponds  with 

V 

a.  Time  of  concentration 

A  most  widely  used  term  is  the  time  of  concentration,  T 
(3,  8,  14,  20,  22)  defined  previously  inChapter  II.  The  Soil  Conser¬ 
vation  Service  relates  T^  to  rise  time  by  the  following  formulas: 

Rise  time  =  0.  5  D  +0.6  T,-, 

e  ^ 
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or 

Rise  time  -  T^,  1  +  0.  fj  T^, 

where  D  is  duration  of  rainfall  excess, 
e 

Clark  (22)  measured  T  as  time  from  end  of  rainfall  excess  to 
point  of  inflexion  on  the  recession  segment  of  the  hydrograph  (figure  9). 
Hence,  his  relationship  between  the  rise  time  and  time  of  concentration 
becomes 

Rise  time  =  D  +  T_ 

e  C 

b.  Time  of  equilibrium,  T  (figure  10) 

For  uniform  rainfall  intensity,  the  time  of  concentration  would 
be  the  time  of  equilibrium  at  which  the  rate  of  runoff  is  equal  to  the 

rate  of  rainfall  supply  (21).  Theoretically,  this  is  the  maximum 

* 

possible  rise  time.  Since  rainfall  excess  does  not  occur  uniformly 
for  any  extended  period  of  time,  equilibrium  flows  are  seldom  en¬ 
countered.  However,  for  small  areas  where  channel  storage  can  be 
neglected  and  the  conditions  approach  two-dimensional  overflow  (unit 
width  concept)  a  time  of  equilibrium  can  result.  Izzard  (23)  showed 
that  time  of  equilibrium  is  a  function  of  channel  length,  rainfall  inten¬ 
sity,  and  slope  and  land  surface.  Morgali  (24)  cited  experimental 
verification  of  the  time  distribution  of  overland  flow  based  on  a  finite 
difference  solution  by  computer  of  the  continuity  and  momentum 
equations  where  boundary  conditions  include  zero  flow  at  upstream 
end,  sub -critical  flow  at  downstream  exit,  and  uniformly  added 


. 
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rainfall  excess  along  the  channel  slope.  The  time  of  equilibrium  is 
represented  as  follows: 


T 

e 


T  0.593  0.605 

c  L _ n _ 

.0. 388  0.  38 

1  ^ 


o 

where  T^=  time  of  equilibrium 
C  =  a  constant 
L  =  length  of  overland  flow 
n  =  Manning  roughness  coefficient 
i  =  intensity  of  rainfall  excess 
S  =  slope  of  surface. 


c.  Lag  time  parameters 

The  best  measure  of  basin  lag  would  be  T  ,  the  time  from 
centroid  of  rainfall  excess  to  centroid  of  the  resulting  runoff  hydro¬ 
graph  (figure  9),  but  usually  this  value  cannot  be  calculated  con¬ 
veniently  (25).  In  view  of  this  difficulty,  various  alternative  concepts 
have  been  introduced.  Snyder  (26)  suggested  that  basin  lag  could  be 
represented  by  Tp ,  the  time  from  centroid  of  rainfall  to  hydrograph 
peak.  Clark  (22)  utilized  T^ ,  time  of  concentration,  as  a  measure 
of  maximum  lag  time  in  the  basin.  Clark's  definition,  in  a  physical 
sense,  is  analogous  to  the  time  of  concentration  of  the  rational  formula 
but  is  independent  of  time  of  rainfall.  Gray  (15)  found  that  lag  time, 

T-q,  is  equal  to  rise  time  T  . 

Jr  K 
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Time  - — * 

Figure  10.  Definition  sketch  of  time  to  equilibrium 


Figure  11.  Definition  sketch  of  storages  defined  by  Laurenson 


Time  - *— 

Figure  12.  Definition  sketch  of  translation  time  defined  by  Bell 
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Laurenson  (27  )  suggested  that  a  watershed  behaves  like  a 

number  of  concentrated  storages  acting  in  series  (figure  11).  Each 

storage  is  expressed  by 

S  =  K  q 
n  n  n 

where 

=  storage  volume 

q  =  outflow, 
n 

He  showed  that  K  has  the  dimensions  of  time  and  is  equal  to  the  time 
between  centroid  of  inflow  hydrograph  and  centroid  of  outflow  hydro¬ 
graph  if  the  storages  are  linear.  When  the  watershed  is  treated  as  a 
whole,  the  time,  ,  between  the  centroid  of  rainfall  excess  and  the 
centroid  of  hydrograph  at  the  outlet  is  given  by: 

K  =  K.  +  K  +  K_  + _ K 

•  12  3  n 

He  suggested  that  this  be  regarded  as  the  "average  storage  delay  time" 
of  the  watershed. 

Bell  (28)  used  two  concentrated  non-linear  storages  to  simulate 
streamflow  behavior  as  shown  in  figure  12.  The  translation  time,  T^, 
is  the  time  between  the  peaks  of  the  inflow  and  outflow  hydrographs  of 
storage  B.  In  practice  it  is  derived  from  the  time,  T  ,  between  the 
end  of  the  main  excess  rainfall  and  the  peak  of  the  outflow  hydrograph 
(figure  9).  It  is  used  for  routing  the  flow  through  storage  B,  i.  e. 
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where 

q  =  outflow 

S  =  volume  of  B  storage 

=  f(q)  =  average  translation  time  from  q^  to  q. 

The  rise  time  is  therefore  equal  to  duration  of  rainfall  excess  plus 
translation  time. 


d.  Base  time 

The  base  time  T  is  the  base  length  of  the  hydrograph  (figure  9) 
In  the  triangular  hydrograph  method  (3)  by  the  Soil  Conservation  Ser¬ 
vice,  Tg  is  related  to  T^  by  the  following  formula: 

tb  =  2-  67  tr 


where 


T  =  0.  5  D  +0.6  T  . 
R  e  C 


3 .  Rise  Time  and  the  Unit  Hydrograph  Concept 

The  unit  hydrograph  was  originated  in  1932  by  Sherman  (7  )  and 
has  since  become  widely  accepted  as  a  method  of  predicting  extreme 
floods,  given  the  volume  of  runoff.  Although  intended  for  streams 
having  at  least  a  few  years  of  stream  flowr  records,  the  method  has 
been  extended  to  ungaged  watersheds  by  means  of  "synthetic  unit 
hydrographs"  and  "dimensionless  hydrograph"  techniques.  At  the 
same  time,  several  researchers  have  extended  and  formalized  the 
theory  mathematically  through  the  use  of  the  "instantaneous  unit 
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hydrograph".  These  developments  have  been  summarized  by  Linsley 
et  al. ,  (25  ) . 

All  of  these  techniques  depend  on  the  basic  unit  hydrograph 
assumption  that  the  discharge  at  any  time  is  proportional  to  the  volume 
of  runoff.  It  follows  that  unit  hydrograph  rise  time  and  other  runoff 
time  parameters  are  constant  for  a  given  watershed  or,  in  other  words, 
are  unaffected  by  the  magnitude  of  the  runoff.  The  principle  of  super¬ 
position  also  follows  the  assumption  that,  if  "excess  rainfall"  or 
runoff  supply  is  divided  into  segments,  the  unit  hydrograph  method 
can  be  applied  to  each  segment  separately  and  the  resulting  hydro¬ 
graphs  can  be  added  together.  Therefore,  the  rise  time  of  the  sum¬ 
mation-curve  varies  only  with  the  storm  duration.  Since  these  are 
characteristics  of  a  linear  system,  the  unit  hydrograph  is,  therefore, 

based  on  assumptions  of  linearity. 

In  recent  years,  these  assumptions  have  been  questioned. 

Knowing  that  the  velocities  of  flowing  water  and  wave  celerities  in  a 
natural  channel  usually  increase  with  discharge,  one  intuitively  expects 
rise  time  and  other  time  parameters  to  be  affected  somewhat  by 
discharge.  Studies  with  laboratory  models  by  Amorocho  (29)  and 
Pabst  (30)  indicate  that  non-linearity  does  exist.  An  analysis  of  field 
data  by  Laurenson  (27)  showed  that  lag  time  varies  with  discharge, 
indicating  non-linearity.  Using  a  mathematical  watershed  model 
Machmeier  (19)  showed  that  rise  time  and  other  time  parameters 
presently  defined  should  not  be  considered  constant  for  a  given 


watershed. 
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4.  Variability  of  Rise  Time  Within  Watersheds 

Some  studies  imply  that  rise  time  within  watersheds  is  approx¬ 
imately  constant  but  the  previous  analysis  shows  that  this  is  difficult 
to  justify.  Unit  hydrograph  rise  times,  concentration  times,  equil- 
brium  times,  lag  times  and  translation  times  are  all  approximately 
constant,  but  the  rise  time  depends  on  the  duration  of  rainfall  excess 
which  varies  from  storm  to  storm  according  to  the  type  of  rainfall, 
antecedent  moisture  and  similar  factors.  Kibler's  study,  discussed 
in  the  literature  review,  showed  a  direct  relationship  between  rise 
time  and  rainfall  duration. 

Rise  time  is  probably  also  affected  by  the  areal  distribution 
of  rainfall.  If  the  storm  is  centered  near  the  outlet,  the  time  of  rise 
should  be  shorter  than  if  the  storm  is  centered  near  the  headwaters. 

Since  there  are  no  previous  studies  of  the  range  of  rise  times 
that  might  be  expected  within  any  particular  watershed,  it  appears 
desirable  that  this  matter  be  examined  in  further  detail. 

5 .  Variability  of  Rise  Time  Between  Watersheds 

Watershed  characteristics  are  sources  of  variability  of  rise 
time  between  watersheds  as  indicated  by  studies  of  basin  lag,  time 
of  concentration  and  the  other  approximately  constant  factors.  Var¬ 
iability  in  rise  time  between  watersheds  can  also  be  studied  directly 
by  selecting  representative  rise  times  such  as  average  or  median. 


This  will  be  attempted  in  the  next  chapters. 


I 
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6.  Conclusions  from  the  Analysis 

The  variability  of  rise  times  within  watersheds  needs  to  be 
analysed  with  a  large  amount  of  data.  If  it  can  be  shown  that  a  general 
distribution  of  values  applies  to  most  watersheds,  this  would  provide 
a  basis  for  the  following  additional  studies: 

1.  Selection  of  a  parameter  or  representative  value  for  com¬ 
paring  rise  times  between  watersheds. 

2.  Derivation  of  a  formula  for  estimating  the  representative 
rise  time  from  watershed  characteristics. 

3.  Evaluation  of  other  methods  of  estimating  rise  time  and 
similar  characteristic  times. 

4.  Possible  use  of  rise  time  for  estimating  design  floods. 
These  will  be  attempted  in  the  remainder  of  the  investigation. 


' 
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CHAPTER  IV 

PRELIMINARY  EXAMINATION  OF  DATA 

1 .  The  Data  Available 

This  study  is  based  on  hydrologic  data  assembled  under  the 
auspices  of  the  "Small  Watershed  Program"  at  Colorado  State  Uni¬ 
versity.  All  watersheds  with  five  or  more  events  were  selected,  thus 
giving  a  total  of  47  watersheds  and  407  flood  events.  The  distribution 
of  watersheds  according  to  their  geographic  locations  by  state  is  shown 
in  table  1.  It  is  noted  that  most  of  the  watersheds  are  from  12  western 
and  mid-western  states  representing  diversified  regions  of  the  drier 
parts  of  the  United  States. 

A  total  of  22  watershed  characteristics  for  each  watershed  was 
taken  out  of  reference  (48).  The  method  of  computing  these 
watershed  characteristics  as  defined  by  Laurenson  et  al. ,  (31  )  is 
shown  in  Appendix  A. 

Six  major  rainfall  and  runoff  parameters  are  recorded  in  table 
B-2  of  Appendix  B.  These  are: 

(a)  antecedent  precipitation  index,  API, 


F* 


(b)  flood  peak,  q, 

(c)  flood  frequency,  P 


. 
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TABLE  1.  DISTRIBUTION  OF  WATERSHEDS  AND  FLOOD 

EVENTS  BY  STATES 


Number  of  Number  of 

Number  States  Watersheds  Events 


1 

Arizona 

8 

49 

2 

California 

6 

33 

3 

Colorado 

1 

5 

4 

Illinois 

1 

5 

5 

Iowa 

1 

5 

6 

Mississippi 

6 

30 

7 

Nebraska 

3 

17 

8 

New  Mexico 

1 

7 

9 

Ohio 

7 

44 

10 

Oklahoma 

2 

14 

11 

Pennsylvania 

1 

140 

12 

Texas 

8 

47 

13 

Wisconsin 

2 

1 1 

TOTAL 

47 

407 
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(d)  total  precipitation,  P  , 

s 

(e)  storm  duration,  D, 

(f)  rainfall  excess  duration,  D  . 

e 

2.  Definitions  of  Rise  Time 

Two  definitions  were  selected  (figures  13  and  14): 

1.  T  ^  is  the  time  from  lowest  discharge  to  peak  of 
the  runoff  hydrograph. 

2.  T  ^  is  the  time  from  discharge  0.  05  (q  -  q^ )  +  q.,  to 
peak  q^  of  the  runoff  hydrograph. 

T  ^  is  usually  implied  by  the  term  "rise  time"  but  sometimes 
this  value  is  a  doubtful  representation  of  the  main  period  of  rise  e.  g. , 
when  there  is  a  long  hydrograph  segment  of  approximately  constant  or 
slightly  increasing  discharge.  In  such  cases,  it  is  also  difficult  to 
decide  when  the  rise  commences. 

T  is  intended  to  avoid  the  above  difficulties.  It  is  essentially 
a  measure  of  the  significant  rise  time.  Both  values  have  been  deter¬ 
mined  from  the  runoff  data  and  are  recorded  in  table  B-2  of  Appendix 
B. 

3.  Correlation  Between  the  Two  Measures  of  Rise  Time 

Figure  15  shows  the  plots  of  T  against  T  For  measuring 

K  1  K  2 

the  degree  of  association  of  the  two  rise  times,  the  correlation  coef¬ 
ficient  was  computed.  It  is  equal  to  0.  90. 


' 

I 
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Figure  14.  Rise  time  definitions  sketch  no.  2 


hours 


30 
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Figure  15.  Correlation  between  T 


and  T 


31 


4.  Correlations  of  Rise  Time  with  Rainfall  and  Runoff  Parameters 

A  preliminary  analysis  of  the  data  was  made  with  simple 

graphical  correlations  of  rise  time  with  rainfall  and  runoff  parameters. 

For  each  watershed,  T  and  T  were  plotted  against  several 

K  1  R  3 

types  of  API,  q,  P  ,  and  D.  No  significant  relationships  were  obtained 

s 

(figure  16).  Most  of  the  watersheds  have  only  five  events,  hence,  most 
of  the  graphs  have  only  five  points. 

With  an  attempt  to  combine  all  events  together,  the  rise  time 
was  converted  to  the  ratio  of  the  actual  rise  time  to  the  mean  rise 
time  of  each  watershed. 

The  ratios  T  /T  and  T  /T  were  plotted  against  API. 

R 1  R 1  RZ  R  Z 

The  drainage  basin  area  was  used  as  a  basis  to  group  the  watersheds 
together  in  the  plotting.  No  significant  relationship  was  obtained 
(figure  17  ) . 

5.  Conclusions  from  the  Preliminary  Study 

T  and  T  are  highly  correlated  with  each  other.  Therefore, 

R  1  R  Z 

T„  ,  will  be  the  only  measure  of  rise  time  considered  in  the  next  cha- 
R 1  J 

pter. 

Due  to  a  lack  of  sufficient  storm  data  pertaining  to  individual 
watersheds,  the  variability  of  rise  times  within  a  watershed  cannot  be 
formulated  in  complete  detail. 

In  an  effort  to  overcome  this  difficulty  and  thereby  obtain  a  gen¬ 
eral  picture  of  rise  time  variability,  it  was  decided  that  further  attempts 
should  be  made  to  combine  the  data  from  different  watersheds. 
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Figure  16.  Correlation  of  rise  time  with  API.  q  ,  P  and  D 
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CHAPTER  V 

THE  CONCEPT  OF  RELATIVE  RISE  TIME 
AND  ITS  DISTRIBUTION 

The  foregoing  chapter  showed  a  need  for  further  attempts  to 
combine  data  of  the  preliminary  study  to  get  a  better  picture  of  varia¬ 
bility  of  rise  times  within  watersheds.  Before  this  can  be  done,  the 
data  has  to  be  transformed  so  that  all  watershed  values  are  comparable 
and  may  be  considered  as  samples  of  a  single  population  of  values. 

1 .  Average  Rise  Time 

First  attempt  at  doing  this  was  to  estimate  the  average  rise 
time  for  all  watersheds  with  five  or  more  events  by  calculating  the 
arithmetic  mean  and  express  actual  rise  times  as  ratios  to  the 
average,  thus  putting  them  in  dimensionless  form. 

It  was  also  considered  that  the  arithmetic  average  rise  time 
would  be  a  good  representative  value  of  rise  time  to  enable  analyses 
of  variability  between  watersheds. 

Much  of  the  data  was  analysed  according  to  these  ideas,  but  it 
soon  became  apparent  that  the  arithmetic  average  rise  time  could  not 
be  accurately  estimated  from  the  small  samples  of  some  watersheds 
because  it  was  unduly  affected  by  extreme  values.  An  attempt  was 
made  to  reduce  this  effect  by  initially  eliminating  rise  time  ratios 
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less  than  0.  5  and  greater  than  2.  0  and  recalculating  the  means.  The 
ratios  of  values  were  then  recalculated  using  these  revised  estimates 
of  the  average.  This  problem  was  later  referred  to  Dr.  Siddiqui  (32), 
and  he  recommended  that  the  median  be  used  rather  than  the  arithmetic 
mean.  In  small  samples  the  median  is  less  affected  by  occasional 
extreme  values. 

2.  Median  Rise  Time 

Median  rise  times  were  calculated  for  all  watersheds  as  shown 
in  table  B-l  of  Appendix  B.  This  parameter  was  adopted  as  the  repre¬ 
sentative  value  of  rise  time  for  each  watershed  and  also  as  a  base 
value  for  converting  the  actual  rise  times  to  dimensionless  relative 
rise  times. 

3 .  Relative  Rise  Times 

"Relative  Rise  Time"  is  defined  here  as  the  ratio  of  actual  rise 
time  to  median  rise  time.  The  relative  rise  times  were  calculated  for 
all  events  as  in  tables  B-2  and  B-3  of  Appendix  B,  and  combined  to 
give  frequency  distributions  for: 

(a)  all  watersheds, 

(b)  a  single  watershed  (Shaver  Creek), 

(c)  watersheds  smaller  than  one  square  mile,  and 

(d)  watersheds  larger  than  one  square  mile. 

These  are  plotted  in  figure  18  which  indicates  similar  distri¬ 


butions  with  positive  skew  for  each. 


- 
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Figure  18.  Histograms  for  frequency  distribution  of  relative  rise  times 
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Figure  19.  Cumulative  probability  distribution  of  relative  rise  times 

(separate  plottings  of  each  sample  ) 
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The  cumulative  frequency  values  were  found  to  be  closely 
linear  when  plotted  on  logarithmic  probability  paper,  indicating  that 
a  logarithmic  transformation  is  suitable  for  normalizing  the  distri¬ 
butions  (figure  19). 

To  test  the  significance  of  differences  between  the  distributions, 
a  Chi-square  (x2)  test  was  applied  with  an  adopted  level  of  significance 
of  five  per  cent.  This  is  described  in  section  4  and  implies  the 
following  hypothesis: 

Hypothesis  -  The  frequency  distribution  obtained  from  the 
combined  records  of  all  watersheds  expresses  the  variability  of  rise 
times  for  any  single  watershed. 

4.  Chi-square  Test  of  Hypothesis 

The  Chi-square  test  is  summarized  as  follows:  The  total  range 

of  sample  observations  is  divided  into  k  class  intervals,  each  having 

the  observed  class  frequency  O.  and  corresponding  expected  class 

1 

probability  E.  (j  =  i,  2,  .  .  . ,  k).  The  measure  of  total  discrepancy 
J 

between  observations  and  expected  values,  X  ,  becomes 
k  lO.  -  E.  1  2 


This  statistic  is  distributed  asymptotically  as  Chi-square  with  (k-1  ) 
degrees  of  freedom  (33). 


' 
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a.  First  test  of  hypothesis:  single  watershed 

The  distribution  of  relative  rise  times  of  a  single  watershed, 
Shaver  Creek,  was  tested  against  the  distribution  of  relative  rise 
times  obtained  from  the  combined  records  of  all  47  watersheds.  The 
computations  are  shown  in  table  B-5  of  Appendix  B. 

Calculated  X2  =  12.  48  . 

For  ten  degrees  of  freedom:  X2  __  =  18.  3  and  X2  =  3.  94. 
&  0.  95  0.  05 

Therefore,  the  difference  is  not  significant. 

b.  Second  test  of  hypothesis:  smaller  than  one  square  mile 
watersheds 

The  46  watersheds,  excluding  Shaver  Creek,  were  separated 
into  two  groups  of  23  watersheds  by  the  size  of  area  for  further  testing 
of  the  hypothesis.  In  this  section,  the  first  group,  smaller  than  one 
square  mile  watersheds,  is  tested.  The  computations  are  shown  in 
table  B-6  of  Appendix  B. 

Calculated  X2  =  13.  75  . 

For  20  degrees  of  freedom:  X2  =  31.4  and  X2  =  10.9. 

6  0.95  0.05 

Therefore,  the  difference  is  not  significant. 

c.  Third  test  of  hypothesis:  larger  than  one  square  mile 
watersheds 

The  second  group,  larger  than  one  square  mile  watersheds, 
is  tested  here.  The  computations  are  shown  in  table  B-7  of  Appendix 
B. 

Calculated  X2  =  21.  53  . 
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For  18  degrees  of  freedom:  X2.  =  28.  9  and  X2  =  9.  39. 

0. 95  0.  05 


Therefore,  the  difference  is  not  significant. 


5.  Theoretical  Distribution  of  Relative  Rise  Times  in  a  Watershed 

Since  a  logarithmic  transformation  has  been  found  suitable  for 
normalizing  the  distribution,  it  is  advantageous  to  derive  its  specific 
mathematical  form. 

The  expression  for  the  log-normal  density  function  with  two 
parameters  is  of  the  form  (33): 


In  x  -  |jl  ]2 


p(x)  = 


1 


e 


o  2 
2  (J 


n 


n 


where 


N 


lnx 


n 


i  -  1  i 
N 


and 


a 


n 


lnx.  -  u  l2 
i  n  I 
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Using  the  graphical  method,  figure  20,  X^  and  were  obtained: 


X  «  \i  =  In  0.  85  =  -  0.  16 
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Figure  20.  Cumulative  probability  distribution  of  relative  rise  times 

(combined  plottings  of  the  four  samples) 
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S  ~  (j  =  In  0.  50  =  -  0.  GO 
n  n 


which  gives  the  following  equation: 


-1.74  R,r  (' 


jin  H,p+  0.  1 6 J 2 
0.  96 


This  is  the  theoretical  distribution  of  relative  rise  times  for 
all  the  watersheds  studied. 


6.  Conclusions  and  Discussion  of  the  Concept 

Rise  time  within  a  watershed,  far  from  being  constant,  varies 
from  about  0.  3  to  3  times  the  median  or  more. 

A  single  distribution  of  relative  rise  times  applies  to  all 
watersheds  investigated.  This  has  the  following  mathematical  form: 


1 

-1.74Rt  6 


jin  Rt+  0.  1 6 J 2 
0.  96 


where  is  the  relative  rise  time. 

From  the  hydrologic  point  of  view,  the  rise  time  or  time  to 
peak  of  the  hydrograph  of  runoff  from  a  drainage  basin  is  influenced 
by  two  major  groups  of  factors: 

(a)  factors  that  are  relatively  constant  for  a  particular  water¬ 
shed  such  as  area,  slope  and  drainage  density; 

(b)  factors  that  vary  with  time  for  a  particular  watershed  such 
as  precipitation  intensity,  precipitation  duration  and  antecedent  moisture. 


\ 
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The  distribution  ol'  relative  rise  times  reflects  variations  due 
to  group  (b).  These  factors  exhibit  seasonal  changes  in  accordance 
with  the  climatic  environment  and  also  random  differences  from  storm 
to  storm.  However,  their  effects  on  the  rise  time  appear  to  be  con- 
sistant  for  a  wide  range  of  conditions  when  expressed  in  the  form  of 
the  distribution  of  relative  rise  times. 

There  may  be  systematic  differences  in  the  distribution  due 
to  geographical  location,  climatic  factors  and  watershed  characteristics 
but  this  study  has  not  been  able  to  detect  any  such  differences  with  the 
available  sample  of  data. 

The  median  rise  time  represents  the  effects  of  the  relatively 
constant  factors  in  group  (a)  and  will  be  studied  in  the  next  chapter. 


' 
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CHAPTER  VI 

PREDICTION  OF  MEDIAN  RISE  TIME 

Several  methods  have  been  suggested  for  predicting  rise  time 
or  time  to  peak,  as  cited  in  the  literature  review.  The  best  known 
and  most  widely  used  in  American  engineering  practice  is  the  Soil 
Conservation  Service  Method  (9).  So,  it  is  considered  worthwhile 
to  examine  the  prediction  by  this  method  before  deriving  other  for¬ 
mulas  expressing  the  relationship  between  rise  time  and  watershed 
characteristics. 

1 .  Evaluation  of  SCS  Method  for  Computing  Rise  Time 

The  SCS  mehtod  for  computing  rise  time  or  time  to  peak  was 
summarized  in  section  3  of  Chapter  II.  Estimated  time  of  concen¬ 
tration,  T  ,  and  rise  time,  T  for  each  watershed  are  shown  in 
C  R 

table  B-l  of  Appendix  B. 

a.  Time  to  peak  (see  tables  B-8  and  B-9  -  Appendix  B) 

Times  to  peak  were  plotted  against  median  rise  times  (figure 
21  ).  A  regression  line  of  best  fit  by  least  squares  was  obtained  and 
has  an  equation  as  follows: 

T  =  0. 656  +  0. 142  T 
1VI  K 

where  T,,  is  the  observed  median  rise  time  and  T  time  to  peak  by 
M  R  J 

SCS  method.  The  coefficient  of  determination,  R2,  is  0.  088.  The 


.  * 
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Figure  21.  Correlation  between  median  rise  time  and  time  to  peak 

of  SCS  method 


46 


average  error  of  estimate  is  145  per  cent.  This  seems  to  agree  also 
with  the  finding  from  the  tests  done  by  the  Hydrological  Research  Unit 
of  the  University  of  Witwatersrand,  Republic  of  South  Africa  (46). 

b.  Time  of  concentration  (see  tables  B-8  and  B-9  -  Appendix  B) 

Times  of  concentration  were  also  plotted  against  median  rise 
times  (figure  22).  The  average  error  of  estimate  of  median  rise  times 
by  times  of  concentration  by  SCS  method  is  103  per  cent.  It  is  an 
improvement  compared  to  the  estimation  of  median  rise  time  by  time 
to  peak. 

c.  Discussion 

From  the  literature  review  it  was  noted  that  the  SCS  method 
is  based  on  ihe  study  of  Kirpirch  who  used  the  data  collected  by 
Ramser  on  six  agricultural  watersheds  ranging  in  area  from  1.  25  to 
112.0  acres,  near  Jackson,  Tennessee. 

The  improvement  of  estimation  of  the  median  rise  time  by  the 
time  of  concentration  is  understandable  because  Ramser  actually 
measured  the  time  of  concentration  by  the  rise  time  which  is  the  time 
required  for  the  water  in  the  channel  at  the  gauging  station  to  rise 
from  the  low  to  the  maximum  stage  as  recorded  by  the  water-stage 
recorders.  The  time  of  concentration  calculated  by  SCS  method  is, 
therefore,  the  average  rise  time  although  this  does  not  appear  to  be 
realized  by  the  SCS  and  other  users  of  the  method. 
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Figure  22.  Correlation  between  median  rise  time  and  the  time  of  concentration 

of  SCS  method 
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It  is  concluded  that  there  is  room  for  improving  the  prediction 
of  rise  time.  Multiple  regression  analysis  will  be  used  to  derive 
formulas  for  estimating  the  median  rise  time. 

2.  Regression  Analysis 

The  regression  and  correlation  analysis  is  one  of  the  oldest 
statistical  tools  used  in  hydrology.  It  was  first  used  for  filling 
missing  data  and  extending  short  records  at  one  hydrologic  station 
by  relating  the  available  data  at  this  station  with  those  at  adjacent 
stations.  Now  its  application  has  been  broadened  to  cover  the  study 
of  the  relationship  between  two  or  more  hydrologic  variables  and  also 
the  investigation  of  dependence  between  the  successive  values  of  a 
series  of  hydrologic  data  (34). 


a.  Multiple  linear  regression 

The  association  of  three  or  more  variables  can  be  investigated 
by  multiple  regression  and  correlation  analysis  (34  and  36). 

The  multiple  regression  relation  may  be  expressed  in  the  form 


where  x^  ,  x^,  .  .  .  ,  x^  are  m  variables.  This  equation  gives  the 
estimates  of  x^  for  given  values  of  all  other  variables.  If  the  rela¬ 
tionship  is  linear,  the  regression  is  called  multiple  linear  regression 
and  the  association  is  multiple  linear  correlation.  Variables  of  non¬ 
linear  relations  in  hydrology  are  often  transformed  to  linear  relations 


, 
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for  multiple  regression  analysis  because  linear  equations  are  easier 
to  treat  than  non-linear  multiple  relations. 

If  there  are  m  variables  to  correlate,  including  one  dependent 
and  m-1  externally  independent,  the  general  equation  for  multiple 
linear  regression  is 


x  =  B  +  B  x  +  .  .  .  +  B.x.  +  .  .  .  +  B  x 
112  2  11  mm 


where  B^  is  the  intercept  and  B^  is  the  multiple  regression  coefficient 
of  the  dependent  variable  x^  on  the  independent  variable  x.,  with  all 
other  variables  kept  constant.  By  the  method  of  least  squares,  m 
linear  equations  may  be  obtained  enabling  one  to  determine  the  m 


unknown  coefficients  as  follows: 


B  S(Ax  ) 2  +  B  2(Ax  Ax  )+...+  B  2(Ax  Ax  )  =  2(Ax  Ax  ) 
2  z  3  z  3  mZrn  iz 

B  S(Ax  Ax  )  +  B  S(Ax  )  2  +  .  .  .  +  B  2(Ax  Ax  )  =  2  (Ax  Ax  ) 
2  23  3  3  m  3  m  1  3 7 


B_  2(Ax_Ax  )  +  B_  2(Ax.Ax  )  +  .  .  .  +  B  2(Ax  )2  =  2(Ax.Ax  ) 
2  2m  3  3m  m  m  1  m 


B  -  x  -  B  x  -  B  x  -  ...-  B  x 
1  1  2  2  3  3  mm 


A. 

where  x.  =  — —  and  Ax.  =  x.  -  x.  ,  with  i  taken  from  1  to  m,  and  N  is 
l  N  ill 

the  sample  size.  These  equations  enable  the  determination  of  m  para¬ 
meters.  The  computational  work  increases  with  the  increase  of  either 
the  number  m  of  the  variables  or  the  sample  size  N. 
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If  the  relationship  is  a  curvilinear  one  of  the  form 


b 


m 


m 


a  logarithmic  transformation  will  yield  a  linear  equation  of  the  form 


log  x ^  =  log  b  +  b2  log  x2  +  b^  log  x^  +  .  .  .  +  b^  log  x 


Inherent  in  the  application  of  multiple  regression  methods  of 
analysis  to  hydrologic  problems  are  four  assumptions  (37): 


1.  there  are  no  errors  in  the  independent  variables;  errors 
occur  only  in  the  dependent  variable; 

Z.  the  variance  of  the  dependent  variable  does  not  depend  on 
the  values  of  the  independent  variables; 

3.  the  observed  values  of  the  dependent  variable  are  serially 
uncorrelated  random  variables; 

4.  the  independent  variables  are  not  correlated  with  each 
other. 

The  multiple  regression  and  correlation  analysis  is  used  a 


great  deal  at  present  because  the  tedious  numerical  computations  can 
be  executed  very  quickly  by  the  aid  of  electronic  digital  computers. 

b.  Computer  program  of  stepwise  regression  analysis 

A  stepwise  multiple  linear  regression  program  (35),  available 
at  the  Engineering  Research  Center  of  Colorado  State  University,  was 
utilized  to  derive  an  equation  for  predicting  a  median  rise  time  of  a 


watershed.  This  program  was  originated  by  the  Health  Sciences 
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Computing  Facility  at  U.  C.  L.  A.  and  was  modified  to  run  on  CDC 
6600  computer. 

This  program  computes  a  sequence  of  multiple  linear  regres¬ 
sion  equations  in  a  stepwise  manner.  At  each  step  one  variable  is 
added  to  the  regression  equation.  The  variable  added  is  the  one  which 
makes  the  greatest  reduction  in  the  error  sum  of  squares.  Equiva¬ 
lently,  it  is  the  variable  which  has  highest  partial  correlation  with  the 
dependent  variable  partial  on  the  variables  which  have  already  been 
added;  and,  equivalently,  it  is  the  variable  which,  if  it  were  added, 
would  have  the  highest  F  value.  The  control  of  the  variable  entered 
at  each  step  is  by  F-value  test,  the  significance  level  being  selected 
by  the  operator.  The  F-values  for  entry  and  deletion  of  variables 
can  also  be  selected  at  any  desired  level  of  significance. 

3.  Application  of  Regression  Analysis 

a.  General 

It  is  worthwhile  at  this  point  to  note  some  of  the  difficulties  of 
the  multiple  regression  technique  when  applied  to  hydrologic  problems. 
The  assumptions,  cited  previously,  upon  which  the  multiple  regression 
approach  is  based  are  often  violated  in  hydrologic  problems. 

The  first  assumption,  that  no  errors  occur  in  the  independent 
variables,  is  violated  in  hydrology.  In  this  investigation,  none  of  the 
variables  can  be  exact  measurements.  However,  it  is  expected  that 
the  dependent  variable  contains  relatively  more  error  than  any  one 
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of  the  independent  variables.  The  second  assumption  is  frequently 
suspect  in  hydrologic  studies,  and  efforts  should  be  made  to  check 
its  significance  in  applications  of  regression  analysis.  In  many 
hydrologic  studies  assumption  three  is  violated,  but  in  this  investi¬ 
gation  the  values  of  rise  times  are  evidently  uncorrelated  random 
variables.  With  regard  to  the  fourth  assumption,  that  the  independent 
variables  are  not  correlated  with  each  other,  the  stepwise  regression 
procedure  usually  results  in  the  first  few  variables  being  almost 
independent  (38  ) . 

The  multiple  correlation  coefficient,  R ,  is  the  primary 
statistical  parameter  used  to  appraise  the  prediction  equations  of 
this  study.  The  square  of  R,  called  the  multiple  coefficient  of  deter¬ 
mination,  is  a  good  indicator  of  the  utility  of  a  derived  equation.  It 
expresses  the  fraction  of  variance  of  the  dependent  variable  attribu¬ 
table  to  the  variation  of  independent  variables  contained  in  the 
regression  equation. 

b.  Estimating  median  rise  time  by  watershed  characteristics 

In  this  part  of  the  study,  only  watershed  characteristics  were 
used  in  the  stepwise  regression  analysis  of  median  rise  times  by 
computer.  Both  values  of  rise  time  measurements  were  used. 

First,  covariance  and  correlation  matrices  of  all  variables 
listed  in  table  2  were  obtained;  then  the  variables  were  reselected 


according  to  their  degree  of  correlation  with  the  median  rise  time  and 
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TABLE  2.  VARIABLES*  USED  IN  THE  REGRESSION  ANALYSIS 


No. 

Symbols 

Names 

Units 

1 

A  - 

Area 

Square  miles 

2 

C 

Compactness  coefficient 

Dimensionless 

3 

Dd 

Drainage  density 

Miles  per  sq.  mi. 

4 

F 

Form  factor 

Dimensionless 

5 

H 

Total  fall 

Feet 

6 

L 

Length  of  main  stream 

Miles 

7 

L 

c 

Distance  to  centroid  of  area 

Miles 

8 

L 

m 

Dimensionless  mean  travel  distance 

Dimensionless 

9 

L 

s 

Total  length  of  extended  streams 

Miles 

10 

Lt 

Average  travel  distance 

Miles 

11 

P 

Perimeter  of  catchment 

Miles 

12 

Ri 

Overland  slope  (1st  definition) 

Feet  per  mile 

13 

R2 

Overland  slope  (2nd  definition) 

Feet  per  mile 

14 

R3 

Overland  slope  (3rd  definition) 

Feet  per  mile 

15 

R4 

Overland  slope  (4th  definition) 

Feet  per  mile 

16 

R5 

Overland  slope  (5th  definition) 

Feet  per  mile 

17 

R6 

Overland  slope  (6th  definition) 

Dimensionless 

18 

S1 

Stream  slope  (1st  definition) 

Feet  per  mile 

19 

S2 

Stream  slope  (2nd  definition) 

Feet  per  mile 

20 

S3 

Stream  slope  (3rd  definition) 

Feet  per  mile 

21 

*4 

Stream  slope  (4th  definition) 

Feet  per  mile 

22 

T 

Ml 

Median  rise  time  (1st  definition) 

Hours 

23 

T 

M2 

Median  rise  time  (2nd  definition) 

Hours 

24 

W 

Average  width  of  catchment 

Miles 

*  The  definitions  of  these  variables  are  in  Appendix  A. 
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their  degree  of  interdependency  between  themselves.  Since  any 
variable  can  be  forced  into  the  equation  in  the  computer  program  as 
mentioned  earlier,  only  variables  best  correlated  with  the  median 
rise  time  and  practically  independent  of  each  other  were  used  to 
obtain  the  prediction  equation. 

The  variables  ennumerated  in  the  equations  to  follow  are  in 
descending  order  of  significance.  Only  statistically  significant 
variables  are  included.  The  F-test  at  a  ten  per  cent  level  was  used, 
as  described  in  the  previous  section  on  computer  program,  to  deter¬ 
mine  whether  an  independent  variable  is  significant  or  not.  The 
physical  meaning  of  each  contribution  was  considered  because 
physical  comprehension  and  interpretation  are  necessary  for  deciding 
the  limitations  and  suitability  of  the  equations  in  practical  application. 
A  summary  of  the  equations  discussed  in  the  ensuing  paragraphs  is 
presented  in  table  3. 

c.  Equations  for  all  watersheds 

Forty-six  watersheds  and  ZZ  catchment  characteristics  listed 
in  table  Z  were  used.  The  results  showed  the  highest  positive  cor¬ 
relations  of  the  median  rise  times  with  the  stream  slopes  and  overland 
slopes.  The  equations  obtained  are  as  follows: 

T__ .  =  -  1.  Z7  +  0.  15  H  x  10_Z  +  0.  57  R  cx  10"2  I’ll 

Ml  5  L  J 


' 

» 

< 


55 


and 

TM2  =  “°-  77  +  °*°8  H  x  10  +  °-  40  R5  x  10  [2] 

where  H  is  the  total  fall  in  feet  and  R_  is  the  overland  slope  in  feet 

5 

per  mile.  T  ^  and  T^.^  are  the  two  measures  of  median  rise  time 
in  hours.  The  coefficient  of  determination  is  0.  75  and  0.  74  res¬ 
pectively  for  equations  1  and  2. 

The  positive  correlations  of  the  median  rise  times  with  the 
slopes  is  contrary  to  physical  reasoning.  Hydraulic  considerations 
suggest  that  the  median  rise  time  should  tend  to  decrease  with 
increases  in  slope.  Graphical  plottings  of  the  median  rise  times 
against  the  stream  slope  and  the  overland  slopes,  figure  23,  show 
that  the  relationship  is  mainly  determined  by  the  six  California 
watersheds  which  differ  greatly  from  the  other  watersheds  in  the 
following  respects: 

1  )  the  flood-producing  storms  are  of  the  long  duration  storms 
(ten  hours  and  longer)  rather  than  the  local  convective 
thunderstorm  type  (three  hours  and  shorter); 

2)  Extremely  steep  topography. 

The  second  assumption  of  regression  analysis  is,  therefore, 
violated  here  because  the  variance  of  median  rise  times  does  depend 
on  the  values  of  the  independent  variables.  It  is  an  unfortunate 
coincidence  that  these  two  extreme  conditions,  1  and  2,  should 
occur  on  the  same  watersheds  as  there  appears  to  be  no  physical 
relationship  between  them.  It  seems  that  the  long  rise  times  are 
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Figure  23.  Median  rise  times  versus  stream  slopes  and  overland  slopes 
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due  to  (1)  rather  than  (2)  but  no  parameters  expressing  (1)  were 
included  in  the  regression  analysis  and  so  the  correlation  is  with  (2). 

This  idea  may  be  tested  by  examining  the  correlation  when  the 
California  watersheds  are  removed  from  the  analysis,  as  will  be  done 
in  the  next  section. 

These  findings  suggest  some  possible  classification  of  water¬ 
sheds  according  to  the  types  of  flood-producing  storms,  i.e.,  long 
or  short  duration  storms  (39).  Detailed  investigations  of  rise  times 
for  long  duration  storms  are  not  justified  here  because  there  are  only 
six  watersheds  from  California  which  fall  within  this  category. 

d.  Watersheds  subject  to  short  duration  storms  only 

This  grouping  only  eliminates  the  California  watersheds 
leaving  40  watersheds.  The  correlation  coefficient  of  rise  times  with 
slopes  becomes  almost  zero.  The  resulting  equations  are  as  follows: 

T  =  0.  42  -  0.  03  D  .  +  0.  48  C  [3] 

Ml  d  l  J 

and 

T  =  0.  34  -  0.  02  D  +  0.  34  C  [Yl 

M2  d  L  J 

where  is  the  drainage  density  in  miles  per  square  mile  and  C  is 
the  dimensionless  compactness  coefficient.  T  ^  and  T^^  are  the 
two  measures  of  median  rise  times  in  hours.  The  coefficient  of 
determination  is  0.  21  and  0.  24  respectively  for  equations  3  and  4. 


- 
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This  indicates  that  only  about  20  per  cent  of  the  variation  in  specific 
rise  time  has  been  explained. 

Logarithmic  transformation  was  also  used  and  the  equations 
became: 


log  TM1  =  0.  16  -  0.44  log  Dd 

and 

l0gTM2=  °-°4  '  °’44  l0gDd 


5 


where  is  the  drainage  density  in  miles  per  square  mile.  T  ^  and 
^M2  are  tw0  measures  °f  median  rise  times  in  hours.  The  co¬ 
efficient  of  determination  is  still  low  0.  25  and  0.  28  respectively  for 
equations  5  and  6. 

In  both  cases  the  coefficient  of  determination  is  too  low  for 
reliable  estimates  and  further  subdivision  appears  necessary.  Exam¬ 
ination  of  hydrographs  suggests  that  there  are  significant  differences 
between  hydrograph  rises  in  humid  areas  and  arid  areas.  The  water¬ 
sheds  are  grouped  in  two  regions:  humid  and  sub-humid  region,  arid 
and  semi-arid  region  as  specified  in  reference  40. 


e.  Humid  and  sub-humid  watersheds  -  short  duration  storms 
For  the  18  humid  and  sub-humid  watersheds  the  equations 
obtained  are  as  follows: 


0.  81  +  0.  36  L 

c 


' 
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and 

TM2  =  °*  86  +  °-  25  Lc  "  °-  65  F  [8] 

where  L  is  the  distance  from  the  outlet  to  centroid  of  area  in  miles 
c 

and  F  is  the  dimensionless  form  factor  (ratio  of  the  area  to  the  square 
of  the  length  of  main  stream,  A/L2).  T  ^  and  T^^  are  two 
measures  of  median  rise  times  in  hours.  The  coefficient  of  deter¬ 
mination  is  0.  38  and  0.  69  respectively  for  equations  7  and  8. 
Logarithmic  transformation  was  also  used: 

L°gTM,--  0.08  +  0.47  log  Lc  [?>] 


and 

L°g  Tm2  =  -  0.  08  +  0.  51  log  Lc 


where  L  is  the  distance  from  the  outlet  to  the  centroid  of  area  in 
c 

miles.  T,,.  and  T  are  the  two  measures  of  median  rise  times 
Ml  M2 

in  hours.  The  coefficient  of  determination  is  0.  55  and  0.  68  res¬ 
pectively  for  equations  9  and  10. 

Sixty-nine  per  cent  of  the  variation  of  the  median  rise  time 
measurement  by  T^^  can  be  explained  by  equation  8,  but  only  38 
per  cent  can  be  explained  if  the  measurement  was  made  by  T 
This  suggests  that  T  is  a  better  measurement  of  median  rise  time 

lVi.  ^ 

than  T  The  result  may  be  explained  by  the  data  which  shows 

that  83  per  cent  of  the  hydrographs  in  the  humid  and  sub-humid  region 


i 

' 
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have  relatively  long  rising  limb,  but  the  initial  part  of  these  limbs  is 
often  a  small  part  of  the  hydrograph  and  not  closely  related  to  the 
main  rise. 


f.  Arid  and  semi-arid  watersheds  -  short  duration  storms 

The  equations  obtained  for  22  arid  and  semi-arid  watersheds 
are  as  follows: 


T 


Ml 


=  1.33  -  0.86  s|x  10  1  -  0.93  A  x  10 


-2 


and 


1.  05  -  0.  68  S2  x  10 


-1 


0.  62  A  x  10 


where  S  is  the  stream  slope  in  feet  per  mile  and  A  is  the  area  in 

Cm! 

miles.  T  and  T^^  are  the  two  measures  of  median  rise  time  in 
hours.  The  coefficient  of  determination  is  0.  68  and  0.  58  respectively 
for  equations  11  and  12. 

Using  logarithmic  transformation,  the  equations  become: 


Log  T 


Ml 


0.  91  -  0.  60  log  S,  -  0.  14  log  L 

3  s 


and 

Log  TA/ro  =  0.  80  -  0.  60  log  S  -  0.  11  log  L 
M2  3  s 


where  is  the  stream  slope  in  feet  per  mile  and  L^  is  the  total  length 
of  extended  streams  in  miles.  T  and  T  are  the  two  measures  of 

1VI  1  -LVJ.  L-t 
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median  rise  times  in  hours.  The  coefficient  of  determination  is  0.73 
and  0.  69  respectively  for  equations  13  and  14. 

The  above  equations  show  that  rise  time  decreases  as  area 
increases  in  arid  regions.  This  agrees  with  the  findings  of  Keppel 
and  Renard  (16),  as  cited  in  the  literature  review,  on  their  hydro¬ 
graph  analysis  for  ephemeral  streams  in  Walnut  Gulch,  Arizona. 
Keppel  and  Renard  claim  that  this  is  due  to  the  high  transmission  loss 
and  the  translatory  wave  effect. 

For  humid  and  sub-humid  regions,  the  correlation  matrix 
shows  the  opposite  effect,  i.e.,  rise  times  tend  to  increase  with 
increasing  areas  and  these  opposing  effects  are  probably  responsible 
for  the  low  coefficients  of  determination  of  equations  3, 4,  5  and  6  which 
include  all  watersheds. 

The  above  results  show  also  that  T,,,  is  a  better  measurement 

Ml 

of  median  rise  time  than  T  ,  but  the  difference  is  not  as  distinct 

IV  J.  ^ 

as  in  the  humid  and  sub-humid  regions.  This  is  reasonable  because 
all  the  available  hydrographs  of  arid  or  semi-arid  watersheds  have 
short  rising  limbs. 


g.  Correlation  between  the  two  measures  of  median  rise  times 

It  has  been  shown  in  the  preliminary  study  that  T  is  highly 

K 1 

correlated  with  T^^.  T,,.  and  T  _  are  also  highly  correlated  and 

R2  Ml  M2  &  J 

their  relationship  is  as  follows: 


T 


Ml 


M2 


0.  12  +  1.  50  T 


ya 

T&J 
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TABLE  3.  RISE  TIMES  PREDICTION  EQUATIONS 


Region 

Sample 

Size 

Equation 

Number 

EQUATIONS 

R2 

F 

All  Watersheds 

46 

1 

2 

T__,  =  -  1. 27+  0.  15  Hx  10~2+  0.  57  R  x  lo’2 
Ml  5 

T,  =  -  0.  77+  0.  08  Hx  10_2+  0.  40  R  x  10~2 
M2  5 

0.75 

0.  74 

63.  6 

60.  3 

TMl  =  -°-12+1-50TM2 

0.  99 

4068 

All  Watersheds 
Short  Duration  Storms 

40 

3 

4 

5 

6 

T  =  0.  42  -  0.  03  D  +  0.  48  C 

Ml  d 

T  =  0.  34  -  0.  02  D  +  0.  34  C 

M2  d 

log  tmi=  °*16"0-44  log  Dd 

log  T'M2=  0.  04  -  0.  44  log  Dd 

0.  21 

0.  24 

0.  25 

0.  28 

5.  0 

5.  8 

12.  5 

14.  9 

TM1=0-01  +  U1  TM2 

0.83 

188 

Humid  and  Sub-Humid 
Watersheds 

18 

7 

8 

9 

10 

TM1=0-81  +  0-36  Lc 

T_  =  0.  86  +  0.  25  L  -  0.  65  F 

M2  c 

l0g  TM1=  0,06+  °‘47  l0g  Lc 

l0g  TM2=  "  °-  08  +  °*  51  log  Lc 

0.  38 

0.  69 

0.  55 

0.  68 

9.  6 

16.  4 

19.  2 

33.  4 

TM1=0-14+1-27  TM2 

0.  75 

46.  7 

Arid  and  Semi-Arid 
Watersheds 

22 

11 

12 

13 

14 

T  =  1.  33  -  0.  8  6  st  x  10_1  -  0.  93  A  x  10  "2 

Ml  2 

T  =  1.  05  -  0.  68  Six  10_1  -  0.  62  A  x  10  "2 

log  TM1  ”  °-91-°-60  logV°- 14  lQg  Ls 

log  T.._  =  0.  80-  0.  60  log  S-  0.  1 1  log  L 

M2  3  s 

0.  68 

0.  58 

0.  73 

0.  69 

20.  6 

13.  0 

25.  4 

21.  1 

TMf°'06+U2TM2 

0.  88 

153 
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for  all  watersheds. 


T  =  0  01  +  1  31  T 
Ml  •  M2 


for  all  watersheds  except  California  watersheds, 


T  =  0  14  +  1  27  T 
Ml  *  M2 


for  humid  and  sub-humid  watersheds,  and 


TM1=  0.°6+  U2T 


for  arid  and  semi-arid  watersheds.  T__.  and  T„.  are  the  two 

Ml  M2 

measures  of  median  rise  times  in  hours. 

The  coefficients  of  determination  are  .75  to  .  99  as  shown 
in  table  3. 


4.  Conclusions  on  Prediction  of  Median  Rise  Time 

The  following  conclusions  can  be  made  from  this  study  of 
rise  time  by  catchment  characteristics: 

(a) .  There  are  a  number  of  difficulties  in  applying  the 
multiple  regression  method  to  hydrologic  data  and  a  high  coefficient 
of  determination  does  not  necessarily  mean  that  an  equation  is 
suitable  for  practical  applications.  The  equations  derived  from  the 
small  watersheds  data  are  summarized  in  table  3. 

(b) .  Equations  derived  for  all  watersheds  were  determined 
mainly  by  the  six  California  samples  which  are  subject  to  long 
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duration  storms.  These  were  not  suitable  for  other  watersheds 
subject  to  short  duration  storms. 

(c) .  For  watersheds  subject  to  short  duration  storms,  the 
rise  time  tends  to  increase  with  area  in  humid  regions  and  decrease 
with  area  in  arid  regions.  These  opposing  effects  are  probably 
responsible  for  low  coefficients  of  determination  of  equations  3,4,  5  and  6 
which  include  all  watersheds. 

(d) .  Reasonable  estimates  can  be  obtained  with  equations  8 
and  10,  which  are  limited  to  humid  and  sub-humid  watersheds  subject 
to  short  duration  storms,  and  equations  11  and  13,  which  are  limited 
to  arid  and  semi-arid  watersheds  subject  to  short  duration  storms. 

The  results  suggest  that  these  equations  give  better  estimates  than 
the  widely  used  SCS  method,  although  it  has  not  been  possible  to 
thoroughly  check  this  with  independent  data. 

(e) .  T  ^  and  T^^  are  highly  correlated  and  the  regression 
equations  are  summarized  in  table  3.  T  is  a  better  measurement 

1VJL  ^ 

of  rise  times  than  T  ^  in  humid  regions.  However,  in  the  arid 
regions  the  two  measurements  are  almost  the  same. 


' 
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CHAPTER  VII 

APPLICATIONS  OF  MEDIAN  AND 
RELATIVE  RISE  TIMES 

For  estimation  of  the  design  flood  hydrograph  it  is  often 
desirable  to  know  what  rise  time  is  appropriate  for  floods  of  various 
frequencies.  There  could  possibly  be  a  tendency  for  lower  rise  times 
to  coincide  with  larger  floods  and  vice  versa. 

There  are  various  hydrologic  techniques  that  require  estimates 
of  median  rise  time,  critical  duration  and  similar  watershed  charac¬ 
teristics.  Some  of  these  will  be  discussed  in  this  chapter. 

1 .  Correlation  of  Rise  Time  with  Flood  Frequency 

If  there  were  many  observations  from  one  watershed,  it  would 
be  desirable  to  investigate  the  variation  in  time  of  rise  with  magnitude 
of  flood  by  correlating: 

-  rise  time  vs.  flood  peak, 

-  rise  time  vs.  volume  of  discharge, 

-  rise  time  vs.  total  rainfall, 
etc .  .  . 

If  some  relationship  could  be  established  between,  say,  rise 
time  and  flood  peak  for  one  watershed,  it  may  be  of  the  form  shown 


by  curve  A  of  figure  24. 


. 
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b  c  >|D0d  P00| J 


Figure  24.  Correlation  between  rise  time  and  flood  peak 
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For  a  larger  watershed,  the  relationship  would  be  expected 
to  be  different  because  rise  times  tend  to  be  larger.  Curve  B  of 
figure  24  might  be  typical  of  a  larger  watershed. 

With  the  small,  watershed  data  there  are  a  few  observations 
from  many  watersheds  and,  if  one  were  to  lump  them  all  together 
to  try  and  get  a  curve  like  A,  one  would  be  unsuccessful  because  of 
the  variations  due  to  watershed  size,  etc.  .  . 

To  remove  the  above  variations,  relative  rise  time  is  sub¬ 
stituted  for  rise  time  and  peak  frequency  for  flood  peak.  Now  all 
observations  are,  essentially,  measures  of  the  deviation  from  the 
average  condition  and  can  be  added  together  to  give  one  large  homo¬ 
geneous  sample  instead  of  a  number  of  small  heterogeneous  samples. 

Figure  25  (a)  shows  plottings  of  relative  rise  times,  R  , 
against  the  corresponding  flood  frequencies.  It  was  concluded  that 
there  is  no  relationship  between  the  relative  rise  times  and  the  flood 
frequencies.  This  result  was  also  checked  with  the  plottings  of  the 
relative  rise  times  versus  flood  frequencies  of  140  events  from  one 
single  watershed  (Shaver  Creek,  figure  25  (b)). 

Since  there  is  no  correlation  between  relative  rise  times  and 
flood  frequencies,  one  can  use  the  median  rise  time  to  estimate  the 
design  hydrograph  for  any  flood  frequency. 
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Flood  Frequency  in  Terms  of  Probobility,  percentage  Flood  Frequecy  in  Terms  of  Probability,  percentage 
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Figure  25.  Correlation  between  relative  rise  time  and  flood  frequency 
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2.  Use  of  Rise  Time  in  Hydrologic  Techniques 

a.  Rational  formula 

The  rational  formula  is: 

Q  =  CIA 

where  Q  is  the  peak  discharge  in  cubic  feet  per  second,  C  is  a  runoff 
coefficient  depending  on  characteristics  of  the  drainage  basin,  I  is  the 
rainfall  intensity  in  inch  per  hour,  and  A  is  the  drainage  area  in  acres 
(8,  21,  47). 

When  using  the  rational  formula,  one  must  assume  that  the 
maximum  rate  of  flow,  owing  to  a  certain  rainfall  intensity  over  the 
drainage  area,  is  produced  by  that  rainfall  which  is  maintained  for  a 
certain  effective  duration.  Up  to  the  present,  this  duration  is  assumed 
to  be  equal  to  the  period  of  concentration  of  flow  at  the  point  under 
consideration  and  referred  to  as  time  of  concentration.  Ramser,  who 
is  responsible  for  the  empirical  formula  for  the  time  of  concentration 
later  made  known  by  Kirpich,  has  determined  the  time  of  concentration 
by  rise  time  or  time  to  peak  (see  Chapter  II).  The  findings  of  this 
study  suggest  that  rise  time  should  be  used  instead  of  time  of  concen¬ 
tration  for  selecting  the  rainfall  intensity  I  in  the  rational  formula. 

b.  Regional  analysis 

The  analysis  of  hydrologic  data  from  relatively  homogeneous 
areas  to  obtain  general  relationships  between  various  hydrologic 


» 
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factors  for  the  region  covering  the  areas  is  sometimes  referred  to 
as  regional  analysis. 

Where  data  are  available  to  develop  flood -frequency  curves 

for  small  watersheds  on  a  regional  analysis,  the  peaks  for  selected 

frequencies  (usually  2-,  10-,  and  100-year  frequencies)  can  be 

plotted  versus  A/T  where  A  is  the  drainage  area  and  T  is  the 

K  K 

rise  time  or  the  time  from  the  beginning  of  rise  to  the  peak  rate  of 
flow.  The  resulting  plotting  will  show  a  rather  consistent  relation¬ 
ship  if  the  physical  characteristics  of  the  watersheds  are  not  too 
variable  (41  ).  The  most  appropriate  rise  time  for  these  estimates 
is  the  median.  The  U.  S.  Geological  Survey  publishes  circulars 
(on  a  state  basis)  which  gives  the  results  of  regional  analyses  of 
watersheds  in  specific  areas. 

c.  Synthetic  hydrographs 

A  synthetic  hydrograph  is  prepared  using  the  data  from  a 
number  of  watersheds  to  develop  a  dimensionless  unit  hydrograph 
applicable  to  ungaged  watersheds  (3). 

1  )  Dimensionless  hydrographs  -  The  curvilinear  hydrograph 
shown  in  figure  26  is  a  dimensionless  unit  hydrograph  prepared  by 
Victor  Mockus  (42)  from  the  unit  hydrographs  of  a  variety  of  water¬ 
sheds.  Various  relationships  between  the  hydrograph  components 
can  be  determined,  the  most  useful  of  which  is  the  equation  for  the 


peak  rate  of  flow: 


« 


' 


7  1 


Figure  26.  Dimensionless  hydrograph  and  mass  curve 

(after  V.  Mockus,  1957) 
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= 


KAQ 


P  T 


R 


where  is  the  peak  rate  of  flow  in  cubic  feet  per  second,  K  is  a 

constant  dependent  on  the  shape  of  the  hydrograph,  A  is  the  drainage 

area  in  square  miles,  Q  is  the  volume  of  runoff  in  inches  and  T  is 

K 

the  rise  time  which  is  the  time  in  hours  from  start  of  rise  to  peak 
rate.  A  nomogram  for  the  solution  of  this  equation  is  given  in  figure 
27  for  K  equal  to  484  as  suggested  by  Mockus.  Single -peaked  hydro¬ 
graphs  can  be  constructed  when  any  three  of  the  four  variables  in  the 
above  equation  are  known  with  the  use  of  figure  26. 

Again,  the  median  rise  time  would  be  the  most  appropriate 

value  of  T  to  use  for  a  particular  watershed. 

R 

2)  Triangular  hydrographs  -  The  curvilinear  hydrograph  of 
figure  26  can  often  be  replaced  by  an  equivalent  triangular  hydrograph 
which  is  more  easily  constructed  and,  for  routing  through  reservoirs 
or  stream  channels,  gives  results  about  as  accurate  as  those  obtained 
using  the  curvilinear  hydrograph  (41  ).  Figure  28  shows  a  triangular 
hydrograph  based  on  the  curvilinear  hydrograph  of  figure  26  where 


T.  =  T_.  +  T  . 
b  R  r 


Since  T^  =  1.67  T  ,  by  construction,  then 


V  2-67  tr 


Area  in  acres 
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Figure  27.  Peak  rates  when  Q  =  1  inch  (after  V.  Mockus,  1957) 
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Figure  28.  Triangular  hydrographs  (after  SCS-USDA,  1957) 
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T  may  be  the  median  rise  time  or  other  specific  rise  times, 

K 

depending  on  the  purpose  of  the  study, 
d.  Design  hydrograph  shape 

The  shape  of  a  design  hydrograph  is  sometimes  obtained  by 
the  use  of  the  dimensionless  unit  hydrograph  (43,  44).  This  idea  is 
that  the  usual  hydrograph  for  storms  of  short  duration  (45), 

Q  =  f  (t,  storm  intensity,  watershed  properties) 

can  be  written  in  dimensionless  form: 

Q/Qc=  f2  (t/tc)  . 

The  latter  expression  will  apply  to  any  watershed,  Q  being  the  dis¬ 
charge  rate  and  the  characteristic  of  the  storm  intensity  -  usually 
the  peak  rate  of  discharge  following  the  storm:  t  ,  sometimes  called 
the  "time  lag",  is  a  characteristic  of  the  particular  watershed.  Usu¬ 
ally  t  measures  the  breadth  of  the  hydrograph;  it  might,  for  example, 

be  the  rise  time  or  time  to  peak,  T  ,  which  is  the  time  elapsed  be- 

R 

tween  the  incidence  of  the  runoff  and  the  maximum  discharge. 

3.  Conclusions  on  Applications  of  Median  and  Relative  Rise  Times 
Relative  rise  time  is  not  correlated  with  flood  frequency; 
therefore,  the  median  rise  time  is  a  suitable  representative  value 
for  any  frequency  flood.  As  a  time  characteristic  of  a  particular 
watershed,  rise  time  is  very  usefull  in  hydrologic  techniques  of  flood 
peak  estimation  and  design  hydrograph  prediction. 
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CHAPTER  VIII 

SUMMARY  AND  CONCLUSIONS 

1 .  Statement  of  Problem 

In  many  instances,  the  hydrologist  is  requested  to  determine 
the  outflow  hydrograph  from  a  basin.  Hence,  he  is  concerned  not 
only  with  the  maximum  discharge  and  the  total  volume  but  also  the 
runoff  distribution  with  respect  to  time. 

The  time  of  occurrence  of  the  maximum  discharge  of  the  out¬ 
flow  hydrograph  from  small  watersheds  (less  than  40  square  miles) 
is  the  subject  of  this  study  and  is  referred  to  as  hydrograph  rise  time 
or  time  to  peak. 

2.  Method  of  Study 

This  investigation  uses  407  flood  events  of  47  watersheds  in 
13  states.  The  main  objective  is  to  study  the  variability  of  rise  times 
and  to  try  and  relate  this  with  climatic  factors  and  physiographic 
factors.  It  is  shown  that  a  general  distribution  of  values  applied  to 
most  watersheds  will  explain  the  change  of  rise  times  caused  by 
climatic  factors  such  as  antecedent  moisture  conditions  of  the  water¬ 
shed,  storm  patterns,  etc.  The  median  rise  time  is  then  a  time 
characteristic  of  a  particular  watershed  which  can  be  estimated  by 
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the  physiographic  watershed  characteristics.  Before  deriving  pre¬ 
diction  equations  of  rise  times  by  using  stepwise  regression  analysis, 
an  evaluation  of  an  existing  popular  method  (Soil  Conservation  Service) 
is  made. 

3 .  Results 

The  complexity  of  small  watershed  responses  is  due  to  the 
heterogeneity  between  individual  small  watersheds  with  respect  to 
climate,  geometric  characteristics  and  other  basin  factors.  The 
relatively  sparse  field  data  from  many  small  watersheds  must  be 
pooled  and  analysed  as  a  large  sample  to  yield  maximum  information 
on  floods  and  their  causative  factors.  This  approach  has  been  used 
in  this  study  by  converting  rise  times  to  dimensionless  relative  rise 
times  (ratio  of  individual  rise  time  to  the  median  rise  time).  It  is 
shown  that  the  same  log-normal  distribution  may  be  adopted  for 
relative  rise  times  in  all  the  watersheds  studied. 

In  the  stepwise  regression  analysis,  it  has  been  found  neces¬ 
sary  to  distinguish  the  watersheds  subject  to  long  duration  storms 
(ten  hours  and  longer)  and  short  duration  storms  or  thunderstorms 
(three  hours  and  shorter).  Furthermore,  the  heterogeniety  between 
individual  small  watersheds  does  not  permit  a  single  general  equation 
for  estimating  rise  times  within  the  short  duration  group.  It  is,  there¬ 
fore,  advisable  to  substratify  the  watersheds  further  into  humid  regions 
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and  arid  regions.  Two  different  equations  have  been  obtained  for 
these  regions  with  about  7  0  per  cent  of  the  variation  explained. 

4.  Conclusions 

Rise  times  within  a  watershed,  far  from  being  constant,  vary 
from  about  0.  3  to  3  times  the  median  or  more.  A  single  distribution 
of  relative  rise  times  R^,  applicable  to  all  watersheds  studied,  has 
the  following  mathematical  form: 

p(R^)  =  probability  density  of  R^  (figure  20) 

(In  Rt+  0.  1 6 J 2 

p(Rx)  =  -1.74  Rt  6  0,96 

The  median  rise  time  in  hours  for  a  humid  watershed  subject 
to  short  duration  storms  can  be  estimated  reasonably  from  the 
following  equation: 

T  =  0.  86  +  0.  25  L  -  0.  65  F 
M  c 

where  L  is  the  distance  to  centroid  of  area  in  miles  and  F  is  the 
c 

dimensionless  form  factor. 

The  median  rise  time  in  hours  for  an  arid  watershed  can  be 
estimated  reasonably  from  the  following  equation: 

T  =  1.  33  -  0.  86  S2xl0_1  -0.093  Axio"2 

lvi.  £ 

where  S  is  the  stream  slope  in  feet  per  mile  and  A  is  the  area  in 

L—i 


square  miles. 
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The  median  rise  time  can  be  used  as  a  time  characteristic 
in  hydrologic  techniques  of  flood  peak  estimation  and  design  hydro¬ 
graph  prediction. 


I 

' 


80 


REFERENCES 


81 


REFERENCES 


1.  Minshall,  N.  E.  ,  "Predicting  Storm  Runoff  on  Small  Experi¬ 

mental  Watersheds,  "  Proc.  ASCE,  V.  86,  No.  HY  8, 
p.  17-38,  Aug.  1960. 

2.  Reich,  B.  M. ,  "Design  Hydrographs  for  Very  Small  Watersheds 

from  Rainfall,  "  CER62BMR41,  Colorado  State  Univ. ,  Fort 
Collins,  Colo.,  July  1962. 

3.  "Section  4:  Hydrology,  "  National  Engineering  Handbook,  Soil 

Conservation  Service,  U.  S.  Dept,  of  Agric.,  Washington, 
D.C.,  pp.  3.  16-1  to  3.  16-22,  1957. 

4.  Gray,  Don  M. ,  "Derivation  of  Hydrographs  from  Small  Water¬ 

sheds  from  Measurable  Physical  Characteristics,  " 

Research  Bulletin  No.  506,  Iowa  State  Univ. ,  Ames,  Iowa, 
June  1962. 

5.  Kinnison,  H.  B.  and  Colby,  B.R.,  "Flood  Formulas  Based  on 

Drainage-basin  Characteristics,  "  ASCE  Trans. ,  Vol.  110, 
pp.  849-904,  1945. 

6.  Brater,  E.  F.  ,  "The  Unit  Hydrograph  Principle  Applied  to  Small 

Watersheds,"  ASCE  Proc.  Vol.  65,  pp.  1  191  -1215,  1939. 

7.  Sherman,  L.  K.  ,  "Streamflow  from  Rainfall  by  the  Unit-graph 

Method,"  Eng.  News-Record,  Vol.  108,  pp.  501  -505,  1932. 

8.  Lloyd-Davis,  D.  E. ,  "The  Elimination  of  Storm  Water  from 

Sewerage  Systems,  "  Inst,  of  Civ,  Engs.,  London  England, 
Minutes  of  Proceedings,  Vol.  164,  p.  4l}  1906. 

9.  "Section  4:  Hydrology,  "  National  Engineering  Handbook,  SCS 

U.  S.  Dept,  of  Agric.,  SCS.  Washington,  D.  C. ,  Section  4, 
Supplement  A,  pp.  3.  15-1  to  3.  15-7,  1957. 

10.  Potter,  W.D.,  "Peak  Rates  of  Runoff  from  Small  Watersheds, 
Hydraulic  Design  Series  No.  2,  Division  of  Hydraulic 
Research,  Bureau  of  Public  Roads,  Washington,  D.  C. , 

April  1961. 


* 


. 


Y 

4 


82 


11.  Chow,  Ven  T. ,  "Hydrologic  Determination  of  Waterway  Areas 

for  the  Design  of  Drainage  Structures  in  Small  Drainage 
Basins,  "  University  of  Illinois,  Engineering  Experiment 
-  Station  Bul  letin  No.  482,  104  p. ,  1962. 

12.  Reich,  B.  M.  and  II iemsfra,  L.  A.  V. ,  "Tacitly  Maximized  Small 

Watershed  Flood  Estimates,  "  Journal  of  Hydraulics  Divi¬ 
sion,  ASCE,  Vol.  91,  No.  HY  3,  Proc.  Paper  4339,  pp. 
217-245,  May  1965. 

13.  Pilgrim,  D.H.,  "Discussions  and  Communications,  "  The 

Journal  of  the  Institution  of  Engineers,  Australia,  Vol.  34, 
No.  3,  pp.  57-59,  March  1962. 

14.  Ramser,  C.E.,  "Runoff  from  Small  Agricultural  Areas,  "  Jour. 

Agr.  Res.,  Vol.  34,  No.  9,  pp.  797-823,  May  1927. 

15.  Gray,  Don  M. ,  "Synthetic  Unit-hydrograph  for  Small  Water¬ 

sheds,  "  ASCE,  HY  4,  Proc.  2854,  pp.  33-54,  July  1961. 

16.  Keppel,  R.V.  and  Renard,  K.G.,  "Transmission  Losses  in 

Ephemeral  Stream  Beds,  "  ASCE,  Vol.  HY  3,  pp.  59-68, 
May  1962. 

17.  Wu,  I-Pai.,  "Design  Hydrographs  for  Small  Watersheds  in 

Indiana,"  ASCE,  HY  6,  Proc.  3694,  pp.  35-66,  Nov.  1963. 

18.  Kibler,  D.  F. ,  "An  Analysis  of  Runoff  from  an  Experimental 

Watershed,"  Penn.  State  Univ.  Dept,  of  Civ.  Eng.,  M.S. 
Thesis,  Sept.  1965. 

19.  Machmeier,  R.  E.,  "The  Effect  of  Runoff  Supply  Rate  and 

Duration  on  Runoff  Time  Parameters  and  Peak  Outflow 
Rates,"  Ph.  D.  Thesis,  Dept,  of  Agr.  Eng.,  Univ.  of 
Minnesota,  University  Microfilms  Inc. ,  Ann  Arbor, 
Michigan,  1966. 

20.  Kirpich,  Z.P.,  "Time  of  Concentration  of  Small  Agricultural 

Watersheds,"  Civil  Engineering  (N.  Y.  ),  Vol.  10,  No.  6. 
p.  362,  June  1940. 

21.  Chow,  Ven  T. ,  'Runoff,  "  Handbook  of  Applied  Hydrology  by 

Ven  Te  Chow,  Section  14,  p.  14-7,  McGraw -Hill  Book  Co. , 
New  York,  1964. 


' 


83 


22.  Clark,  C.O.,  "Storage  and  the  Unit  Hydrograph,  "  Trans.  -Am. 

Soc.  of  Civil  Engrs. ,  Vol.  110,  pp.  1419-1488,  1945. 

23.  Izzard,  C.F. ,  "Hydraulics  of  Runoff  from  Developed  Surplus,  " 

Proc,  Hwy.  Research  Board,  Vol.  26,  pp.  129-150,  1946. 

24.  Morgali,  James  R.,  "Hydraulic  Behavior  of  Small  Watersheds,  " 

Dept,  of  Civil  Engr.,  Tech.  Report  No.  30,  Stanford  Univ. , 
1963. 

25.  Linsley,  Kohler  and  Paulhus,  "Hydrology  for  Engineers, 

McGraw-Hill  Book  Company,  Inc.,  New  York,  1958. 

26.  Snyder,  F.  F. ,  "Synthetic  Unitgraphs,  Trans.  Am.  Geophys. 

Union,  Vol.  19,  Part  1,  pp.  447-454,  1938. 

27.  Laurenson,  E.M. ,  "A  Catchment  Storage  Model  for  Runoff 

Routing,"  Journal  of  Hydrology,  Vol.  2,  pp.  141  -163,  1964 

28.  Bell,  F.C.,  "Translation  Time,  "  University  of  N.  S.  W.  Water 

Research  Laboratory  Report,  No.  85,  1967. 

29.  Amorocho,  J.,  "Measures  of  Linearity  of  Hydrologic  Systems,  " 

Journal  of  Geophysical  Research,  Vol.  68,  No.  8,  pp. 
2237-2249,  April  15,  1963. 

30.  Pabst,  A.  F. ,  "Hydrograph  Linearity  in  an  Elementary  Channel,  1 

M.  S.  Thesis,  Department  of  Civil  Engineering,  University 
of  Minnesota,  1966. 

31.  Laurenson,  E.M.,  Yevjevich,  V.  M. ,  and  Shulz,  E.  F. ,  "Re¬ 

search  Data  Assembly  for  Small  Watershed  Floods,  " 
CER63EML-EFS-VMY  37  ,  Engineering  Research  Center, 
Colorado  State  University,  Fort  Collins,  Sept.  1963. 

32.  Siddiqui,  M.  M. ,  Professor  of  Mathematics  at  Colorado  State 

University,  Private  communication,  Nov.  1966. 

33.  Markovic,  R.  D. ,  "Probability  Functions  of  Best  Fit  to  Distri¬ 

butions  of  Annual  Precipitation  and  Runoff,  "  Hydrology 
papers.  No.  8,  Colorado  State  University,  Fort  Collins, 
August  1965. 

34.  Yevdjevick,  V.  M. ,  "Regression  and  Correlation  Analysis,  " 

Handbook  of  Applied  Hydrology  by  V.  T.  Chow,  Section 
8-11,  McGraw-Hill  Book  Company,  New  York,  1964. 


r 


. 


q  .3  'uri  jhjb  M  n->  (v  ^  .K  1  oe  ■  - 

V/ 

■ 


84 


35.  Program  BMD02R,  "Stepwise  Regression,  "  Version  of  April 

9,  1964,  Health  Sciences  Computing  Facility,  U.  C.L.A. 

36.  Spiegel,  Murray,  "Theory  and  Problems  of  Statistics,  "  Schaum's 

Outline  Series,  Schaum  Publishing  Company,  New  York,  1961. 

37.  Snyder,  W.M.,  "Some  Possibilities  for  Multivariate  Analysis 

in  Hydrologic  Studies,  "  Journal  of  Geophysical  Research, 

Vol.  167,  No.  2,  pp.  721-731,  Feb.  1962. 

38.  Wallis,  James  R.,  "Multivariate  Statistical  Methods  in  Hydro¬ 

logy  -  A  Comparison  Using  Data  of  Known  Functional 
Relationships,  "  Water  Resources  Research,  Vol.  1,  No.  4, 
pp.  447-461,  4th  quarter,  1965. 

39.  Pilgrim,  D.H.,  "Flood  Producing  Storms  on  Small  Rural 

Catchments,  "  Institution  of  Engineers,  Australia,  Civil 
Engineering  Transactions,  CE  8,  No.  1,  April  1966. 

40.  Reich,  Brian  M. ,  "Design  Floods  for  Small  Watersheds  in  the 

Arid  and  Semi-arid  West,  "  Colorado  State  University, 

Civil  Engineering  Department,  unpublished  progressive 
report  for  period  of  September  10,  1965  through  June  8, 

1966. 

41.  Ogrosky,  Harold  O.  and  Mockus,  Victor,  "Hydrology  of 

Agricultural  Lands,  "  Handbook  of  Applied  Hydrology  by 
Ven  T.  Chow,  Section  21,  p.  21-41  to  45,  McGraw-Hill 
Book  Company,  New  York,  1964. 

42.  Mockus,  Victor,  "Use  of  Storm  and  Watershed  Characteristics 

in  Synthetic  Hydrograph  Analysis  and  Application,  "  U.  S. 
Conservation  Service,  1957. 

43.  Bender,  D.  L.  and  Roberson,  J.A.,  "The  Use  of  a  Dimensionless 

Unit  Hydrograph  to  Derive  Unit  Hydrographs  for  Some 
Pacific  Northwest  Basins,  "  Journal  of  Geophysical  Research, 
66,  521  -527,  1961. 

44.  Bender,  D.  L. ,  "Unit  Hydrograph  Variation  with  Storm  Hydro¬ 

graph  Shape,  "  Ph.  D.  Thesis,  University  of  Wisconsin, 

Aug.  1963. 

45.  Lienhard,  John  H. ,  "A  Statistical  Mechanical  Prediction  of  the 

Dimensionless  Unit  Hydrograph,  "  Journal  of  Geophysical 
Research,  Vol.  69,  No.  24,  pp.  5231  -5238,  December  15, 
1964. 


f 


* 


a ' 


85 


46.  Midgley,  D.C.  and  Schultz,  G.A.,  "Progress  in  Flood  Hydro¬ 

graph  Synthesization,  "  Die  Siviele  Ingenienr  in  Suid-Afrika, 
Feb.  1965. 

47.  Kuichling,  D. ,  "The  Relation  Between  the  Rainfall  and  the 

Discharge  of  Sewers  in  Populous  Districts,  "  Trans. 

ASCE,  Vol.  20,  pp.  1-56,  1889. 

48.  "Research  Data  Assembly  for  Small  Watersheds  -  Part  II", 

Engineering  Research  Center,  Colorado  State  University, 
Fort  Collins,  Colorado  ( in  preparation  for  publication  - 
expected  date  of  release:  September  1967  ). 


N 


86 


APPENDIX  A 

METHOD  OF  COMPUTING  WATERSHED  CHARACTERISTICS 

This  has  been  extracted  from  an  unpublished  report  entitled: 
"Explanation  of  Symbols  and  Method  for  Filling  in  the  Watershed 
Card"  in  the  Research  Data  Assembly  for  Small  Watershed  Floods 
at  Colorado  State  University,  prepared  by  R.  Markovic,  R.  Downer 


and  M.  Siddiqi,  May  1964. 
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LIST  OF  SYMBOLS 


A. 


C 

Dd 

F 

Ah 

h 

H 

K 

1. 

1 

L 

L 

L 


L 


con 


L 


L 


m 


L 


Area 

Contour  interval 
Compactness  coefficient 
Drainage  density 
Form  factor 
Difference  in  Elevation 
Elevation  above  outlet 
Total  fall 

Hourly  depletion  ratio 
Length  of  reach 

Distance  from  grid  intersection  to  outlet 
Length  of  main  stream 

Average  length  of  two  successive  contour 
lines 

Distance  to  centroid  of  area 

Length  of  a  contour  line 

Length  of  the  grid  lines  within  the  catch¬ 
ment 

Dimensionless  mean  travel  distance 
Total  length  of  extended  streams 


Square  miles 
Feet 

Dimensionless 

Miles  per  sq.  mi. 

Dimensionless 

Feet 

Feet 

Feet 

Dimensionless 

Miles 

Miles 

Miles 

Miles 

Miles 

Miles 

Miles 

Dimensionless 

Miles 
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L 

t 

Average  travel  distance 

Miles 

n 

Number  of  points 

N. 

1 

Number  of  intersections  of  the  grid  lines 
with  a  contour  line 

P 

Perimeter  of  catchment 

Miles 

PA 

Average  annual  precipitation 

Inches 

qA 

Areal  average  discharge 

Cfs  per  sq.  mi 

Qa 

Average  discharge 

Cfs 

R1 

Overland  slope 

Feet  per  mile 

R2 

Overland  slope 

Feet  per  mile 

R3 

Overland  slope 

Feet  per  mile 

R4 

Overland  slope 

Feet  per  mile 

R5 

Overland  slope 

Feet  per  mile 

R6 

Overland  slope 

Dimensionless 

S1 

Stream  slope 

Feet  per  mile 

S2 

Stream  slope 

Feet  per  mile 

S3 

Stream  slope 

Feet  per  mile 

S4 

Stream  slope 

Feet  per  mile 

St 

Standard  deviation 

Miles 

Sd 

Dimensionless  standard  deviation 

Dimensionless 

w 

Average  width  of  catchment 

Miles 

X 

Latitude  of  gaging  station 

Degrees 

X 

r* 

Latitude  of  centroid 

Degrees 

«• 
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Y  Longitude  of  gaging  station 

Longitude  of  centroid  of  area 
Z  Elevation  of  gaging  station 

Z  Average  elevation  of  a  reach 


Degrees 

Degrees 

Feet-mean  sea  level 
Feet 


I 

*  . 

* 

. 


90 


CATCHMENT  CHARACTERISTICS 


Scale 

Since  all  dimensions  are  to  be  expressed  in  feet  or  miles,  it  will  be 
advantageous  to  convert  the  map  scale  to  feet  and  miles  and  draw  it 
on  the  margin  of  the  map. 

Area,  A 

The  area  will  usually  be  given,  but  if  not,  it  will  be  necessary  to  de¬ 
lineate  the  catchment  boundary  and  planimeter  the  area.  Record  on 
the  map  in  some  convenient  place  (e.  g.  below  the  name  of  watershed) 
the  area  in  square  inches  and  the  area  in  square  miles  to  three  deci¬ 
mal  places.  If  the  area  is  given,  one  needs  only  record  the  area  in 
square  miles  on  the  map.  For  the  present,  areas  less  than  64  acres 
and  greater  than  25,  600  acres  will  be  omitted. 

Length  of  Main  Stream,  L 

Extend  all  marked  stream  systems  up  to  the  watershed  boundaries  in 
accordance  with  the  contours.  Exceptions  to  this  rule  will  be  streams 
which  appear  to  originate  in  springs  or  swamps.  The  "main  stream" 
is  defined  as  that  stream  draining  the  greatest  area.  Using  the  paper 
strip  method,  mark  off  the  length  of  the  main  stream  and  at  the  same 


time  mark  and  label  points  where  the  main  stream  crosses  a  contour 
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line.  Calculate  the  distances  in  miles  to  two  decimal  places  between 
successive  contours  and  record  on  the  paper  strip.  The  summation 
of  these  distances  is  the  length  of  the  main  stream,  L.  Save  the  strip 
of  paper  for  use  in  plotting  the  main  stream  profile,  and  calculation 
of  slope. 

Total  Length  of  Extended  Streams,  L 

s 

Using  the  paper  strip  method,  measure  the  total  length  of  all  extended 
streams.  It  is  suggested  to  first  measure  all  the  tributaries  on  one 
side  of  the  main  stream,  then  the  other  side.  The  task  can  be  made 
easier  by  summing  distances  from  key  stream  forks  rather  than  mea¬ 
suring  all  distances  to  the  outlet. 

Length  to  Center  of  Area,  L 

The  centroid  of  the  catchment  can  be  found  quickly  and  easily  and  with 
a  fairly  high  degree  of  accuracy  by  centering  over  a  map  of  the  water¬ 
shed  a  clear  plastic  overlay  with  a  system  of  lines  drawn  on  it  at  45° 
angles  to  form  a  star-shaped  design.  L^  is  the  distance  along  the 
main  stream  from  the  outlet,  to  a  point  adjacent  to  the  center  of  area 
projected  to  the  main  stream.  This  distance  can  be  found  most  easily 
by  using  the  paper  strip  originally  used  to  measure  the  length  of  the 

main  stream,  L.  Determine  L  in  miles  to  two  decimal  places. 

c 

Mean  Travel  Distance,  L^ 

The  average  travel  distance  is  determined  by  measuring  the  travel 


distance  to  the  outlet  along  the  stream  system  from  each  intersection 


* 

•* 


\v' 


92 


of  a  square  grid  placed  over  a  map  of  the  catchment.  The  grid  should 

be  oriented  in  a  North-South,  East-West  orientation.  The  grid  should 

be  of  such  a  size  that  between  30  and  50  intersections  fall  within  the 

catchment.  It  is  suggested  to  number  each  of  the  intersections  which 

fall  within  the  catchment  to  aid  in  accounting  for  every  intersection. 

The  distances  should  be  measured  by  the  paper  strip  method  and 

recorded  in  miles  in  tabular  form  so  that  L  can  easily  be  determined 

21  1 

from  -  ,  where  21  is  the  sum  of  individual  travel  distances  and 

n  t 

n  is  the  number  of  intersections  within  the  catchment.  By  noting  and 
recording  on  the  map  the  distances  to  key  points  or  stream  forks,  the 
measurement  of  distances  to  other  points  can  be  greatly  facilitated. 

Standard  Deviation  of  Travel  Distance,  s 

The  standard  deviation  of  the  travel  distance  will  be  computed  as  the 
square  root  of 


This  form  lends  itself  to  easy  computation  since  the  second  term  is 


equal  to  the  average  travel  distance,  ,  squared.  The  first  term. 
21^2  ,  can  be  determined  on  the  desk  calculator  by  using  accumulative 


multiply. 

Perimeter,  P 

The  perimeter  is  the  distance  around  the  catchment  measured  along 


the  watershed  boundary.  Using  the  paper  strip  method,  determine 


. 

' 

«• 


. 


the  perimeter  in  miles  to  two  decimal  places,  starting  at  the  gage 
and  proceeding  around  the  area  and  back  to  the  gage  to  form  a  closed 
circuit. 

Total  Fall,  H 

Using  the  strip  of  paper  with  the  main  stream  marked  off  on  it  as  the 
abscissas,  plot  a  graph  of  distance  vs.  elevation  along  the  main  stream 
on  20  x  20  to  the  inch  graph  paper.  After  the  profile  is  plotted,  extra¬ 
polate  each  of  the  ends  of  the  profile.  Determine  the  minimum  and 
maximum  elevations  of  the  main  stream  from  these  extended  slopes. 
The  total  fall  can  now  be  determined  to  the  nearest  foot.  Record  the 
distances  in  miles  between  successive  contour  lines  on  the  profile. 

Stream  Slope.  S^ 

S^  is  calculated  by  dividing  the  total  fall,  H  ,  by  the  Length  of  the 
Main  Stream  in  miles,  L  .  Calculate  to  the  nearest  foot  per  mile. 

Stream  Slope,  S^ 


221. z.  221. z. 

S  =  - — —  =  - i— i-  ,  feet/  mile 

21  2  L2 

i 

where  1.  =  distance  along  the  main  stream  between  successive  con¬ 
tours.  The  individual  1  can  be  easily  determined  from  the  plotted 
profile  of  the  main  stream  or  from  the  paper  strip  used  to  measure 
the  length  of  the  main  stream,  z.  =  the  average  elevation  above  the 
outlet  for  each  reach  of  length,  1.  . 
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Stream  Slope, 


(V 


3/2 


(Ah) 


1/2 


where  1  =  length  of  a  reach  in  miles  and  Ah  =  the  change  in  elevation 
in  a  reach.  Calculate  S^  to  the  nearest  foot  per  mile. 

Stream  Slope,  S^ 

g  _  elevation  at  85%  of  stream  length-elevation  at  10%  of  stream  length 

4  75%  of  stream  length 

The  elevation  at  85%  and  10%  of  the  stream  length  can  most  easily  be 
found  by  drawing  vertical  lines  on  the  main  stream  profile  at  85%  and 
10%  of  the  length  and  noting  the  elevation  where  these  lines  cut  the 
profile. 

Overland  Slope,  R 

cL 

R  =  — >  feet/ mile 

1  A 

where  c  =  the  contour  interval.  L  =  the  total  length  of  contour 

con 

lines  on  a  map  of  the  catchment,  and  A  =  the  catchment  area.  L 

con 

can  best  be  determined  by  using  the  paper  strip  method.  Label  each 
strip  and  mark  the  length  of  each  contour  line  on  it  for  later  use. 

„  1.  57  cN 

R2=  — 

g 
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where  c  =  the  contour  interval,  N  =  the  number  of  intersections  of 
the  grid  lines  with  contour  line,  and  L  =  the  total  length  of  grid  lines 
within  the  catchment  measured  in  both  the  North-South  and  East-West 
directions.  L  can  be  determined  using  the  paper  strip  method.  N 
can  be  determined  by  starting  at  the  catchment  boundary  and  following 
each  contour  line  and  counting  the  number  of  times  each  crosses  a 
grid  line. 

Overland  Slope, 

.  2(Ah  x  L) 

R3 - A  ' 

where  Ah  =  the  difference  in  elevation  between  contours,  L  =  the 

average  length  of  two  successive  contours,  and  A  =  the  catchment 

area.  L  can  be  calculated  from  L  found  while  calculating  R,  . 

con  1 

Consider  the  outlet  and  extreme  boundary  point  of  the  watershed  as 
having  contour  lines  of  zero  length  when  computing  values  of  L  . 


Overland  Slope, 

R^  is  the  mean  overland  slope.  The  overland  slope  is  calculated  by 
dividing  the  contour  interval  by  the  perpendicular  distance  between 
contours  at  each  intersection  of  the  square  grid. 

2— 

R  =  — —  =  — —  =  - - — r  ,  feet /mile 

4  n  n  mean  1 


, 
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where  c  =  contour  interval,  1  =  the  distance  between  contours  at  each 
grid  intersection  and  n  =  the  number  of  grid  intersections  for  which  1 
was  computed. 

Overland  Slope,  R 

-  5 

Rj.  is  the  median  of  the  overland  slopes  computed  above.  This  can  be 
found  by  arranging  the  distances  between  contours  in  descending  order. 
The  median  is  the  distance  which  evenly  splits  this  descending  array, 
with  one  half  of  the  distances  above  it  and  one  half  of  the  distances 
below  it.  When  the  median  falls  between  two  distances,  the  average 
of  those  distances  is  the  median. 

Relief  Ratio,  R„ 

-  6 

The  longest  dimension  of  the  basin  is  the  longest  distance  between  two 

parallel  lines  touching  the  boundaries  of  the  catchment. 

Total  Fall,  H  ..  .  1 

R„  =  - - - — - - : - — - : - —  ,  dimensionless 

6  Longest  dimension  of  the  basin 

Calculate  to  4  decimal  places. 

Drainage  Density,  D^ 

The  Drainage  Density,  D^  ,  equals  the  total  length  of  Extended 
Streams,  L^  ,  divided  by  the  Catchment  Area,  A. 

L 

D  =  r  ,  miles /square  miles 

Q 

Average  Width  of  Catchment,  W 

The  Average  Width  of  the  Catchment,  W  ,  equals  the  Catchment  Area, 

A  ,  divided  by  the  Length  of  the  Main  Stream,  L . 


' 

■ .  ■ 


97 


W  =  A  /L  ,  miles 


Form  Factor,  F 

The  Form  Factor,  F  ,  equals  the  Average  Width  of  the  Catchment, 
W  ,  divided  by  the  Length  of  the  Main  Stream,  L  . 

F=  ~r~=  “j~2"  ,  dimensionless 

JLj  J_j 

Compactness  Coefficient,  C 

The  Compactness  Coefficient,  C  ,  equals  0.  28  times  the  Perimeter, 
P  ,  divided  by  the  square  root  of  the  Catchment  Area,  A  . 

.  0.Z8P 

C  =  — == —  ,  dimensionless 

yT 


Dimensionless  Mean  Travel  Distance,  L 
-  m 


The  Dimensionless  Mean  Travel,  L  ,  equals  the  Mean  Travel  Dis¬ 
co 


tance,  L^  ,  divided  by  the  square  root  of  the  Catchment  Area,  A  . 


L 

L  =  ~r=-  ,  dimensionless 

m  Va 


Dimensionless  Standard  Deviation,  S 

-  d 

The  Dimensionless  Standard  Deviation,  S^  ,  equals  the  Standard 
Deviation,  S^  ,  divided  by  the  square  root  of  the  Catchment  Area,  A 

st 

S  ,  =  -frr-  ,  dimensionless 
d  a/A 


Average  Hourly  Depletion  Ratio,  K 

Plot  the  recession  of  all  available  flood  hydrographs  on  semi- 
logarithmic  paper,  discharge  (log-scale)  vs.  time  (arithmetic  scale) 
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Fit  the  straight  line  by  eye  to  each  recession  line  and  read  the  dis¬ 
charges  at  the  end  (Q^)  an^  the  beginning  (Q^  )  of  time  increment  (At) 
equals  one  hour  in  upper  part  of  each  recession.  Compute  the  ratio 

Q  /Q  for  each  of  these  and  then  average  them.  Use  the  semi- 
^  1 

logarithmic  paper  2  cycles  x  70  divisions  (No.  358-61;  K+  2). 


, 
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APPENDIX  B 

DATA  AND  COMPUTATIONS 


TABLE  B- 1 .  TIME  PARAMETERS  IN  HOURS 
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TABLE  B-3.  DATA  FOR  SHAVER  CREEK 

3.  75  sq.  mi. 

Central  Pennsylvania 


Storm 

No. 

Date 

Rise  Time 
Tri,  hrs. 

h=1r1 

Tri 

Peak 

(cfs) 

Frequency 

% 

1 

3/10/43 

13.  0 

1.  15 

42.  90 

0.  24 

2 

5/  6  /44 

15.  0 

1.  33 

87.  00 

0. 

48 

3 

7/12/44 

7.  0 

0.  62 

1.70 

NA 

6 

4/17/45 

13.  0 

1.  15 

11. 50 

0. 

10 

7 

5/17/45 

7.  0 

0.  62 

103. 00 

0. 

57 

8 

5/28/45 

7.  0 

0.  62 

22.  40 

0. 

14 

9 

7/  5 / 45 

4.  0 

0.  36 

3.  13 

NA 

10 

9/17/45 

23.  0 

2.  04 

17. 00 

0. 

12 

11 

11/21/45 

10.  0 

0.89 

33.  00 

0. 

18 

12 

11/28/45 

14.  0 

1. 24 

75.  00 

0. 

42 

13 

3/  9  / 46 

8.  0 

0.71 

35.  40 

0. 

20 

14 

3/26/46 

7.  0 

0.  62 

15.40 

0. 

11 

15 

5/26/46 

17.  0 

1.51 

130.  00 

0. 

69 

16 

6/  2/46 

12.  0 

1. 06 

72.  50 

0. 

41 

17 

7/  1/46 

9.  0 

0.80 

23.  00 

0. 

14 

18 

9/29/46 

7.  0 

0.  62 

4.  96 

NA 

19 

10/12/46 

9.  0 

0.80 

13.  50 

0. 

10 

20 

10/20/46 

13.  0 

1.  15 

7.  20 

NA 

21 

4/16/47 

12.  0 

1. 06 

22.40 

0. 

14 

23 

7/27/47 

11.0 

0.  98 

8.70 

NA 

24 

11/  8/47 

10.  0 

0.89 

10.  50 

0. 

11 

25 

11/11/47 

9.  0 

0.80 

19.  10 

0. 

13 

26 

4/24/48 

8.  0 

0.  71 

20.  00 

0. 

13 

27 

6/12/48 

9.  0 

0.80 

5.  70 

NA 

28 

7/13/48 

10.  0 

0.89 

7.  30 

29 

8/  5/48 

12.  0 

1. 06 

3.  25 

30 

7/18/48 

7.  0 

0.  62 

7.  50 

f 

31 

11/  6/48 

8.  0 

0.  71 

19.  00 

0. 

12 

32 

11/19/48 

24.  0 

2.  13 

16.  70 

0. 

11 

33 

12/  6 / 48 

10.  0 

0.89 

19.  60 

0. 

12 

34 

1/  5/49 

20.  0 

1. 78 

33.  70 

0. 

19 

35 

2/19/49 

13.  0 

1.  15 

18.60 

0. 

12 

36 

3/23/49 

15.  0 

1.  33 

7.75 

0. 

16 

37 

4/13/49 

13.  0 

1.  15 

29. 00 

0. 

17 

38 

5/  2/49 

11.  0 

0.  98 

12.  40 

0. 

10 
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TABLE  B-3.  DATA  FOR  SHAVER  CREEK  -  Continued 


Storm 

No. 

Date 

Rise  Time 
TR1,  hrs. 

R=> 

Tri 

Peak 

(cfs) 

Frequency 

% 

39 

5/22/49 

7.  0 

0.62 

19.40 

0.  12 

40 

7/13/49 

11.  0 

0.  98 

11.85 

0.  10 

41 

7/17/49 

7.  0 

0.  62 

28.70 

0.  16 

42 

8/29/49 

12.  0 

1. 06 

4.80 

NA 

43 

10/30/49 

13.  0 

1.  15 

7. 45 

NA 

44 

5/18/50 

19.  0 

1. 69 

12.  30 

0.  11 

45 

6/  4/50 

15.  0 

1.33 

15.  50 

0.  11 

46 

7/  5/50 

4.  0 

0.  35 

10.  20 

0.  10 

47 

7/10/50 

7.  0 

0.  62 

6.  05 

NA 

48 

8/  2/50 

4.  0 

0.  35 

7.  30 

NA 

50 

10/28/50 

10.  0 

0.89 

16.30 

0.  11 

51 

11/26/50 

25.  0 

2.  22 

397. 00 

0.  99 

52 

12/  4/50 

12.  0 

1. 06 

71. 30 

0.  40 

53 

3/30/51 

8.  0 

0.71 

152. 00 

0.  76 

54 

4/22/51 

15.  0 

1.  33 

32.  10 

0.  18 

55 

4/30/51 

9.  0 

0.80 

25.  00 

0.  15 

56 

9/  6  / 5 1 

7.  0 

0.  62 

3.  75 

NA 

57 

11/14/51 

6.  0 

0.  53 

1. 50 

1 

58 

11/16/51 

8.  0 

0.  71 

1.96 

1 

59 

12/  5 / 5 1 

13.  0 

1.  15 

15.  10 

0.  11 

60 

2/  3 / 52 

18.  0 

1.60 

26.00 

0.  15 

61 

3/11/52 

14.  0 

1. 24 

150. 00 

0.  76 

62 

4/  5/52 

19.  0 

1. 69 

29.60 

0.  17 

63 

4/25/52 

13.  0 

1.  15 

19.  10 

0.  12 

64 

5/12/52 

12.  0 

1. 06 

15.  20 

0.  11 

65 

7/18/52 

3.  0 

0.  27 

4.  70 

NA 

66 

8/18/52 

4.  0 

0.  35 

3.  50 

NA 

67 

12/10/52 

16.  0 

1.42 

56.  30 

0.  31 

68 

3/  4/53 

10.  0 

0.89 

13.  40 

0.  11 

69 

3/15/53 

8.  0 

0.71 

15.  30 

0.  11 

70 

3/24/53 

14.  0 

1. 24 

155. 00 

0.  77 

71 

4/  7  /  53 

13.  0 

1.  15 

15.  50 

0.  11 

72 

4/18/53 

15.  0 

1.33 

15.60 

0.  11 

73 

11/22/53 

6.  0 

0.  53 

4.  96 

NA 

74 

12/  7/53 

17.  0 

1.51 

5.  59 

NA 

75 

2/20/54 

13.  0 

1.  15 

15.  10 

0.  11 

76 

3/  1/54 

10.  0 

0.89 

170. 00 

0.  13 

77 

6/  4/54 

7.  0 

0.  62 

19.  60 

NA 
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TABLE  B-3.  DATA  FOR  SHAVER  CREEK  -  Continued 


Storm 

No. 

Date 

Rise  Time 
TR1,  hrs. 

R  .1*1 
Tri 

Peak 

(cfs) 

Frequency 

% 

78 

7/  4/54 

7.  0 

0.  62 

1.62 

NA 

79 

10/15/54 

13.  0 

1.  15 

155. 00 

0.  77 

80 

12/14/54 

13.  0 

1.  15 

22.  40 

0.  15 

81 

4/21/54 

8.  0 

0.71 

22.  90 

0.  15 

82 

6/  8 / 55 

17.  0 

1.51 

14.  00 

0.  11 

83 

7/  1/55 

5.  0 

0.  44 

2.  30 

NA 

85 

8/13/55 

17.  0 

1.51 

43.  50 

0.  24 

86 

8/18/55 

14.  0 

1.24 

5.90 

NA 

88 

9/23/55 

12.  0 

1. 06 

2.  20 

NA 

89 

10/14/55 

13.  0 

1.  15 

90.  00 

0.  50 

90 

11/16/55 

9.  0 

0.80 

34.  10 

0.  20 

91 

2/25/56 

13.  0 

1.  15 

30.80 

0.  18 

92 

5/  2/56 

11.  o 

0.98 

18.  30 

0.  12 

93 

5/  6/56 

6.  0 

0.  53 

31.40 

0.  18 

94 

6/23/56 

7.  0 

0.  62 

43.  70 

0.  24 

95 

7/  2/56 

13.  0 

1.15 

31. 90 

0.  14 

98 

9/  6/56 

9.0 

0.80 

20.  00 

0.  13 

99 

10/  4 / 56 

13.  0 

1.  15 

4.46 

NA 

100 

10/22/56 

14.  0 

1.24 

49.  60 

0.  17 

101 

11/  2/56 

14.  0 

1.24 

18.70 

0.  12 

102 

11/22/56 

10.  0 

0.89 

8.  20 

NA 

103 

12/14/56 

18.  0 

1.60 

24.  00 

0.  15 

104 

4/  2/57 

13.  0 

1.  15 

28.40 

0.  16 

105 

4/26/57 

4.  0 

0.  35 

36.  10 

0.  20 

106 

10/24/57 

11. 0 

0.  98 

2.  00 

NA 

107 

12/26/57 

6.  0 

0.  53 

32.  70 

0.  19 

ip-*- 

o 

CD 

6/  1 / 58 

10.  0 

0.89 

5.82 

NA 

110 

7/22/58 

21.  0 

1.86 

3.  00 

NA 

111 

8/  3/58 

6.  0 

0.  53 

10.  05 

0.  09 

112 

9/21/58 

17.  0 

1.51 

4.  65 

NA 

113 

3/  6 / 59 

11.  0 

0.  98 

12.  50 

0.  10 

114 

4/10/59 

8.  0 

0.71 

15.  00 

0.  11 

115 

8/  5 / 59 

5.  0 

0.  44 

3.  10 

NA 

116 

10/  1  / 59 

13.  0 

1.  15 

4.  90 

NA 

118 

11/  6 / 59 

9.  0 

0.80 

29.  00 

0.  17 

119 

11/27/59 

12.  0 

1. 06 

80.  00 

0.  44 

120 

12/12/59 

21.  0 

1.86 

38.00 

0.  21 

121 

1/  2/60 

12.  0 

1. 06 

24.  00 

0.  15 

A 

• 

' 
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TABLE  B-3.  DATA  FOR  SHAVER  CREEK  -  Continued 


Storm 

Date 

Rise  Time 

R  -  Tr1 

Tri 

Peak 

Frequency 

No. 

TR1,  hrs. 

(cfs) 

% 

122 

11/  1/60 

13.  0 

1.  15 

5.40 

NA 

123 

11/29/60 

17.  0 

1.51 

4.  90 

NA 

124 

2/26/61 

11.0 

0.  98 

209. 00 

0.  88 

125 

4/16/61 

8.0 

0.  71 

37. 00 

0.  21 

126 

4/22/61 

8.  0 

0.  71 

33.  00 

0.  19 

127 

5/16/61 

8.  0 

0.  71 

16.  70 

0.  12 

128 

7/13/61 

2.0 

0.  18 

3.  50 

NA 

129 

7/15/61 

4.  0 

0.  35 

11. 00 

0.  10 

132 

11/16/61 

16.  0 

1.42 

6.  00 

NA 

133 

11/23/61 

16.  0 

1.42 

15.00 

0.  11 

134 

4/  7 / 62 

23.  0 

2.04 

73.00 

0.  41 

135 

8/  9 / 62 

3.  0 

0.  27 

3.  70 

NA 

137 

4/18/63 

12.  0 

1.06 

8.  50 

NA 

138 

4/29/63 

14.  0 

1.24 

13.  60 

0.  11 

139 

5/18/63 

14.  0 

1.24 

44.  00 

0.  20 

140 

7/  2/63 

4.  0 

0.  35 

2.  70 

NA 

141 

8/  4/63 

4.  0 

0.  35 

4.  30 

1 

142 

11/  6 / 63 

13.0 

1.15 

7.  50 

143 

11/23/63 

11.  o 

0.98 

2.  10 

I 

144 

11/29/63 

15.  0 

1.33 

17.60 

0.  12 

145 

12/  8 / 63 

15.  0 

1.33 

4.  10 

NA 

146 

3/10/64 

13.  0 

1.  15 

105. 00 

0.  58 

147 

8/  3 / 64 

9.  0 

0.80 

10.  00 

0.  09 

148 

5/21/46 

9.  0 

0.80 

35.70 

0.  15 

150 

12/12/49 

9.  0 

0.80 

42.  00 

0.  23 

151 

11/20/50 

15.  0 

1.  33 

20.  10 

0.  13 

152 

12/  7 / 50 

11.  0 

0.  98 

54.  00 

0.  30 

153 

5/  8  /  54 

11.  o 

0.98 

23.  00 

0.  15 

154 

4/17/54 

12.  0 

1. 06 

25.  00 

0.  15 

2  =  1576. 0 

Total  Number  of  Storms 

=  140 

T-  1576.0 

tR1-  140  -11.26 

hrs. 

NA  - 

Not  available 

» 

> 

■ 

' 
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TABLE  B-4.  DATA  FOR  COMPUTATION  OF 
RELATIVE  RISE  TIMES  DISTRIBUTION 

All  Watersheds 


Order 

Number 

Class 

Interval 

Observed 

Frequency 

(0) 

Cumulative 

Frequency 

Cumulative 

Probability 

1 

0.  11-0.  30 

8 

8 

0.  020 

2 

0.  31-0.  50 

27 

35 

0.  086 

3 

0.  51-0.  70 

57 

92 

0.  226 

4 

0.  71-0.  90 

78 

168 

0.  413 

5 

0.  91-1.  10 

82 

252 

0.  619 

6 

1. 11-1. 30 

58 

310 

0.  762 

7 

1.31-1.50 

31 

341 

0.838 

8 

1.  51-1.  70 

27 

368 

0.  904 

9 

1. 71-1. 90 

12 

380 

0.  934 

10 

1. 91-2.  10 

5 

385 

0.  946 

11 

2.  11-2.  30 

4 

389 

0.  956 

12 

2.  31-2.  50 

4 

393 

0.  966 

13 

2.  51-2.  70 

3 

396 

0.  973 

14 

2.  71-2.  90 

2 

398 

0.  978 

15 

2.  91-3.  10 

2 

400 

0.  983 

16 

3.  11-3. 30 

0 

17 

3. 31-3. 50 

0 

18 

3.  51-3.  70 

2 

402 

0.  988 

19 

3. 71-3. 90 

3 

405 

0.  995 

20 

5.  51-5.  70 

1 

406 

0.998 

21 

7. 31-7. 50 

1 

407 

1. 000 

407 

.  1  ! 


• 

TABLE  B-5.  DATA  FOR  COMPUTATION  OF 
RELATIVE  RISE  TIMES  DISTRIBUTION 
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cr 
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aJ 

X! 
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TABLE  B-6.  DATA  FOR  COMPUTATION  OF 
RELATIVE  RISE  TIMES  DISTRIBUTION 
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TABLE  B-7.  DATA  FOR  COMPUTATION 
RELATIVE  RISE  TIMES  DISTRIBUTION 
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TABLE  B-8.  DATA  FOR  EVALUATION  OF  SCS  METHOD 
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TABLE  B-9.  DATA  FOR  EVALUATION  OF  SCS  METHOD 
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r2  =  f  NSxy  -  (  S  x)(  S  y)  ]  2 

[  N  S  x  2  -  (2  x)2]  [  NS  y2  -  (2  y)2] 

[  39  (  62.  87)  -  (  62.  32)  (34.  46)  )2 

[  39  (154.  43)  -  (62.  32  )  2J  [  39  (43.  086)  -  (34.  46)  2]  =  °- 0879 
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BLM  Library 
Denver  Federal  Center 
Bldg.  50,  OC-521 
P.O.  Box  25047 
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